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PREDICTION  OF  SUPERSONIC  STORE 
SEPARATION  CHARACTERISTICS 

Volume  I.-  Theoretical  Methods  and 
Comparisons  with  Experiment 


1.  INTRODUCTION 

This  report  presents  the  results  of  a combined  theoretical/ experi- 
mental research  program  aimed  at  extending  the  subsonic  six-degree-of- 
freedom  store  separation  trajectory  prediction  method  of  references  1 
and  2 to  supersonic  speeds.  The  accompanying  computer  program  is 
described  in  Volume  II  of  this  report,  reference  3.  Data  from  the 
experimental  part  of  the  program  are  available  in  reference  4 . 

The  theoretical  work  reported  herein  is  the  result  of  the  initial 
effort  towards  building  up  a supersonic  store  separation  trajectory 
method  and  computer  program  which  will  eventually  have  the  same  capabili- 
ties as  the  subsonic  method  and  program  described  in  references  1 and  2. 
With  this  view,  the  input  and  output  of  the  supersonic  computer  program 
have  been  designed  to  resemble  the  existing  subsonic  program  to  the  maxi- 
mum possible  extent.  Following  the  subsonic  approach,  the  present  and 
future  work  are  designed  to  proceed  from  simple  configurations  to  more 
complicated  ones  in  flow  modeling  and  store  load  calculation  methods. 

The  present  computer  program  is  designed  to  handle  a configuration  con- 
sisting of  an  axisymmetric  fuselage  and  a wing  without  dihedral  with  one 
pylon  per  half -wing  or  one  pylon  under  the  fuselage.  Both  the  wing  and 
pylon  may  have  thickness  and  the  wing-fuselage  junction  may  be  off  the 
midplane.  This  configuration  can  be  viewed  as  a basic  one  to  which  other 
components  can  be  attached  and  included  in  the  flow  model  systematically 
as  was  done  during  the  advanced  phases  of  the  subsonic  work.  As  such, 
the  present  program  i3  the  logical  starting  point  for  one  which  will 
ultimately  predict  store  separation  trajectories  from  aircraft  with 
features  such  as  a noncircular  fuselage,  air  inlets,  and  MER  and  TER 
store  groupings. 

The  experimental  program  carried  out  in  conjunction  with  the  theo- 
retical work  was  directed  towards  providing  data  to  aid  in  the  development 
of  the  flow  models  for  the  wing,  pylon,  and  fuselage  components.  In 
order  to  isolate  pylon  effects,  data  were  taken  with  and  without  the 
pylon.  The  data  consist  of  flow  field  surveys,  store  pressure 
distributions,  and  store  forces  and  moments. 


The  next  section  of  this  report  describes  in  considerable  detail 
the  flow  models  based  on  linear,  potential  flow  theory  for  the  components 
of  the  supersonic  aircraft  and  the  store.  A special  section  is  devoted 
to  nonlinear  aspects  of  the  flow  field;  namely,  effects  associated  with 
blunt  wing  leading  edges  and  shock  waves  that  do  not  coincide  with  Mach 
lines.  Following  this  the  store  force  distribution  and  total  force  and 
moment  calculation  methods  are  summarized  and  certain  significant  aspects 
are  discussed.  In  connection  with  the  loading  calculation  methods  and 
the  trajectory  calculation  method,  reference  is  made  to  the  applicable 
subsonic  work  of  references  1,  2,  and  5. 

Finally,  comparisons  are  made  between  theory  and  results  obtained 
during  the  wind-tunnel  test  program  in  order  to  assess  the  accuracy  of 
the  theoretical  methods  in  predicting  flow  fields  about  a wing-fuselage 
configuration  with  and  without  the  pylon.  Comparisons  are  shown  for 
three  Mach  numbers  and  two  angles  of  attack.  For  seme  of  these  conditions, 
comparisons  are  presented  between  t.he  measured  and  predicted  store  load 
distributions  and  forces  and  moments.  In  addition,  flow  field  comparisons 
are  made  using  a different  wind-tunnel  model.  Problem  areas  requiring 
further  attention  are  indicated  and  recommendations  are  made  to  improve 
the  present  method. 

2.  GENERAL  APPROACH 

It  is  possible  to  determine  the  forces  and  moments  acting  on  an 
external  store  in  the  presence  of  a parent  aircraft  in  the  following 
manner.  Consider  first  the  entire  airframe  excluding  the  store  in  question 
The  nonuniform  flow  field,  induced  by  the  flow  model  representing  the  air- 
frame, is  computed  in  the  region  occupied  by  the  store.  If  the  store  is 
now  placed  in  this  nonuniform  flow  field,  then  the  calculated  loading 
acting  on  it  includes  primary  interference.  The  primary  interference 
can  be  viewed  as  the  first  term  in  an  iterative  procedure.  Higher-order 
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jected  to  the  primary  interference  flow  field,  its  effects  on  the  airframe 
can  be  determined.  These  effects,  in  terms  of  velocities  normal  to  the 
wing,  pylon,  and  fuselage,  must  then  be  cancelled  by  an  additional  set 
of  singularities  distributed  over  these  components.  A new  nonuniform 
flow  field  can  now  be  computed  in  the  region  occupied  by  the  store,  and 
the  loading  acting  on  it  can  be  recalculated.  At  this  stage,  the  result 
is  associated  with  the  second  term  in  the  iteration  procedure.  The 


difference  between  the  results  of  the  second  and  the  first  loading 
calculation  is  due  to  additional  interference.  This  process  can  be 
continued  until  the  effects  of  additional  interference  become  negligible. 

For  the  purpose  of  expediting  the  procedure  mentioned  above,  the 
approach  taken  here  for  supersonic  flow,  which  was  used  during  the  earlier 
subsonic  work  (refs.  1,  2,  and  5},  involves  a modification  to  the  first 
iteration  after  which  the  calculation  is  stopped.  The  modification  con- 
sists of  including  in  the  first  iteration  the  effects  of  the  store  in  the 
boundary  condition  applied  to  the  airframe  flow  model.  In  general,  the 
flow  field  under  the  wing  tends  to  align  itself  with  the  wing.  Since  the 
store  diameters  are  usually  small  compared  to  the  local  wing  chord,  the 
induced  sidewash  and  downwash  at  points  on  the  wing  due  to  crossflow  over 
the  store  are  very  small.  Volume  effects  can  be  important  in  connection 
with  the  shock  associated  with  the  store  nose.  At  this  point,  the  calcu- 
lated store  loads  can  be  looked  upon  as  the  result  of  one  and  one-half 
iterations.  The  success  of  this  method  depends  on  the  degree  to  which 
the  primary  interference,  accounted  for  in  this  way,  encompasses  the 
significant  effects.  Methods  for  calculating  additional  interference 
have  been  developed  in  reference  5 for  the  subsonic  case.  It  is  shown 
that  the  effect  due  to  this  interference  on  store  loading  is  3mall, 
provided  that  important  interference  causing  components  such  as  the 
pylon  are  included  as  part  of  the  airframe. 

For  the  supersonic  case,  store  volume  effects  are  expected  to  be 
significant  if  the  store  nose  shock  intersects  the  wing  and  reflects 
back  onto  the  store.  In  this  situation,  the  nonuniform  flow  field  in 
which  the  store  is  actually  submerged  may  be  heavily  influenced  by  the 
store  nose  shock  reflected  off  the  wing.  It  is  for  this  reason  that  the 
store  volume  effects  should  be  included  in  the  first  iteration.  An 
important  advantage  realized  by  this  method  is  that  the  need  for  resolving 
the  complete  mutual  interference  problem  between  the  separated  store  and 
airframe  throughout  the  trajectory  is  eliminated,  thereby  saving  large 
amounts  of  computer  time.  If  the  store  nose  shock  influences  the  airframe 
appreciably,  the  method  has  a disadvantage  in  that  the  airframe  loading 
must  be  recomputed  as  the  store  drops  away.  At  present,  the  computer 
program  does  not  recalculate  the  aircraft  loadings  during  a store  trajec- 
tory, but  means  for  simply  overcoming  this  difficulty  are  subsequently 
ri commended . 
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The  flow  models  used  to  represent  the  aircraft  components  and  store 
in  supersonic  flow  will  be  discussed  and  the  geometric  parameters  of  the 
various  components  and  the  flight  conditions  described.  The  fuselage 
is  flow  modeled  first  and  its  effects  on  the  wing-pylon  combination 
determined.  Certain  terms  in  the  flow  tangency  boundary  conditions  are 
explained  in  relation  to  the  mutual  interference  between  the  aircraft 
components.  Shocks  associated  with  the  wing  leading  edge  receive  special 
attention  because  of  their  influence  on  the  interference  flow  field. 
Methods  for  determining  the  store  forces,  moments,  and  trajectories  are 
summarized. 

3 . FLOW  MODELS 

In  general,  the  flow  models  for  the  components  to  be  described  are 
based  on  linear  theory  for  steady,  inviscid  flow.  Because  the  flow 
regime  considered  here  is  supersonic,  the  associated  potential  flow  is 
governed  by  the  wave  equation.  In  a rectangular  coordinate  system 
(x,y,z)  with  the  x-axis  in  the  direction  of  the  free-stream  vector,  the 
potential  0 is  given  by 

(M^  ~ 1)  ~ m 0 (1) 

dx2  dy£  dz2 

This  equation  is  valid  if  the  perturbation  velocities  u,v,w  given 
below  in  terms  of  the  potential  are  small. 


Because  of  the  linearity  of  the  potential  0 in  equation  (1) , it  is 
possible  to  make  use  of  superposition.  This  principle  allows  for  the 
representation  of  the  aircraft  components  by  distributions  of  potentials 
or  singularities.  Flow  conditions  at  a given  field  point  are  given  by 
the  summed  influences  of  all  the  singularities.  In  particular,  the 
influence  of  each  singularity  is  included  in  the  flow  tangency  condition 
as  will  be  shown  in  a later  discussion  concerned  with  the  flow  tangency 


condition.  Before  describing  the  flow  models  in  detail,  the  geometric 
characteristics  of  the  airframe  and  store  are  listed  next. 

3.1  Description  of  Aircraft  Components  and  Flight  conditions 

The  configuration  of  interest  consists  of  a wing  attached  to  a 
fuselage.  A pylon  can  be  positioned  under  the  wing  or  under  the  fuselage 
centerline. 

The  geometric  parameters  of  the  aircraft  components  described  below 
are  accounted  for  in  the  computer  program.  Flight  conditions  cover  the 
supersonic  speed  regime  and  are  specified  below  together  with  the 
expected  range  of  validity  of  the  present  method. 

3.1.1  Fuselage 

The  fuselage  is  axisymmetric  and  has  a pointed  nose.  A fuselage 
meridian  is  described  by  a set  of  polynomials.  The  afterbody  need  not 
be  pointed. 

3.1.2  Wing 

The  wing  may  have  twist  and  camber  and  the  leading  and  trailing 
edges  can  have  breaks  in  sweep.  Dihedral  angle  and  incidence  relative 
to  the  fuselage  centerline  are  not  accounted  for.  vertical  location  of 
the  wing  attachment  on  the  fuselage  is  not  limited  to  the  midplane. 
Effects  of  wing  thickness  are  included  in  the  program.  The  wing  leading 
edge  can  be  blunt  or  sharp. 

3.1.3  Pylon 

Leading  and  trailing  edges  may  be  swept  without  breaks  in  sweep. 

The  length  of  the  pylon  root  chord  need  not  equal  the  local  wing  chord. 
Thickness  is  accounted  for  and  the  leading  edge  can  be  blunt  or  sharp. 

The  pylon  can  be  located  under  the  wing  or  fuselage  but  is  aligned  with 
the  free  stream. 

3.1.4  Store 

The  general  shape  is  a body  of  revolution  with  a pointed  nose.  A 
body  meridian  is  described  by  a set  of  polynomials.  The  store  can  have 
one  set  of  planar  or  cruciform  stabilizing  fins.  The  present  program 
will  only  handle  a single  store.  Initial  velocity,  pitch  angle  and 
position,  and  the  initial  rolling,  pitching,  yawing  rates  can  be 
prescribed.  .'he  store  may  be  powered  with  a specified  thrust-time 
history.  Store  mass  and  moments  of  inertia  are  constant  and  the  center 
of  gravity  may  be  located  off  the  longitudinal  axis. 
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3.1.5  Flight  conditions 


The  flight  path  is  straight,  angle  of  attack  is  constant,  and  angle 
of  sideslip  is  zero.  The  Mach  number  should  lie  between  1.2  and  3.0. 

The  angle  of  attack  should  be  in  the  0°-10°  range. 

The  flow  model  for  each  of  the  components  will  be  described  next. 

It  will  be  shown  how  the  mutual  interference  between  the  wing,  pylon, 
and  fuselage  is  accounted  for. 

3.2  Flow  Model  for  the  Fuselage  and  Store  Body 

The  potential  flow  method  used  to  represent  an  axisymraetric  body  at 
incidence  in  supersonic  flow  is  discussed  in  this  section.  The  basic 
theory  is  given  in  references  6,  7,  and  8.  Fundamentally,  an  axisym- 
metric  body  in  supersonic  flow  can  be  represented  by  a distribution  of 
line  sources  or  sinks  and  doublets  on  the  body  centerline  to  account 
for  volume  and  angle  of  attack  effects,  respectively.  The  strengths 
of  these  singularities  are  determined  from  the  flow  tangency  conditions 
applied  at  points  on  the  body  surface.  Consider  a body  of  revolution 
and  the  cylindrical  coordinate  system  used  in  the  solution  as  shown 
in  the  following  sketch.  The  angle  between  body  centerline  and  free- 


stream  direction,  a^,  is  assumed  small.  Velocity  component  is 

tangent  to  the  body  contour  in  the  crossflow  plane  and  velocity  component 
vB  is  positive  in  the  radially  outward  direction.  It  is,  therefore, 
necessary  only  to  consider  the  flow  tangency  in  the  meridian  plane 
because  w£  is  already  tangent  to  the  body  surface.  In  any  meridian 
plane  the  body  contour  is  given  by 

rB  “ (3) 

The  flow  tangency  boundary  condition  at  the  body  surface  is 
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The  governing  potential, equation  (l),is  now  expressed  in  the  cylindrical 
coordinate  system  shown  in  the  sketch  above  as 
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Because  of  the  linear  character  of  this  equation,  the  complete  solution 
9 can  be  written  as  the  sum  of  two  potentials. 
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which  satisfies  equation  (5).  These  solutions  are  called  the  axial  and 
crossflow  potentials,  respectively.  The  perturbation  velocities  are 
related  to  the  two  potentials  by 
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Equation  (4)  can  therefore  be  split  into  two  parts  for  bodies  of 
revolution  as  follows. 
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“ c— — ; crossflow  (8b) 
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It  is  noted  that  equation  (8a)  is  the  boundary  condition  for  the 
axial  flow  problem  with  free-stream  velocity,  V^.  Similarly,  equa- 
tion (8b)  involves  only  crossflow  velocities  and  is  the  boundary  condition 
for  the  crossflow  problem  with  crossflow  velocity,  V^a^.  The  sum  of 
equations  (8a)  and  (8b)  yields  the  complete  flow  tangency  condition , equa- 
tion (4) . Thus,  the  first  term  on  the  right-hand  side  of  equation  (6) 
together  with  equation  (8a)  constitute  the  axial  flow  problem.  The 
second  term  on  the  right-hand  side  of  equation  (6)  together  with  equa- 
tion (8b)  formulate  the  crossflow  problem.  These  two  problems  are 
uncoupled  and  can  therefore  be  solved  separately.  Adding  the  two 
solutions  will  satisfy  equations  (4)  and  (5)  exactly. 

The  choice  of  singularity  types  or  potentials  for  the  axial  and 
crossflow  problems  is  described  in  some  detail  in  Appendix  I.  This 
appendix  also  shows  the  incorporation  of  the  boundary  conditions  (8a) 
and  (8b)  in  the  determination  of  the  singularity  strengths  for  the  two 
problems.  A distribution  of  supersonic  cone  solutions  along  the  body 
axis  will  be  used  to  solve  the  axial  flow  problem.  The  crossflow  problem 
is  solved  by  means  of  a distribution  of  supersonic  line  doublets  along 
the  body  axis.  Collectively,  these  singularity  distributions  form  the 
flow  model  for  a body  of  revolution  in  supersonic  flow  based  on  linear, 
potential  flow  theory.  Presently,  the  nonlinear  effects  of  the  shock 
associated  with  the  fuselage  nose  or  store  nose  are  not  accounted  for 
in  the  computer  program.  Instead,  the  first  effects  of  the  fuselage  or 
store  when  traversing  along  a line  parallel  to  the  centerline  are  based 
on  linear  theory  and  are  determined  by  the  Mach  cone  from  the  nose. 

3.3  Plow  Modal  for  Wing-Pylon  Combination  and  Interference  Shell 
on  the  Fuselage 

This  section  describes  the  representation  of  the  wing-pylon  in 
terms  of  elements  from  linear,  potential  flow  theory.  Also,  the  methods 
used  to  account  for  mutual  interference  between  the  wing-pylon  and  the 
fuselage  are  discussed.  Specifically,  constant  u-velocity  panels  are 
distributed  over  the  -wing,  pylon,  and  part  of  the  fuselage.  These  panels 


90c 

SrT  + afC°S  9 

B 


K 


8 


were  originally  called  constant  pressure  panels  by  woodward  in  refer- 
ence 9.  In  addition,  constant  source  panels  distributed  over  the  wing 
and  pylon  account  for  thickness  effects.  The  constant  u-velocity  panels 
laid  out  on  the  fuselage  form  the  interference  shell  and  serve  primarily 
to  cancel  interference  effects  from  the  wing-pylon  combination. 

3.3.1  Geometrical  layout 

The  left  wing  and  pylon  and  the  left  half  of  the  fuselage  are 
divided  up  into  area  panels  of  trapezoidal  shape.  The  right  half  of  the 
airframe  is  accounted  for  by  employing  symmetry  properties  discussed 
later.  Figure  1 shows  a swept  wing-swept  pylon  combination  attached  to 
a fuselage  above  its  midplane.  The  rectangular  coordinate  systems  used 
m the  layout  part  of  the  computer  program  and  associated  with  the  fuse- 
lage (Xg, yB»zB)  and  the  wing  (x^, yw,zw)  are  shown.  The  origin  of  the 
former  is  at  the  fuselage  nose.  The  wing  coordinate  system  has  its 
origin  in  the  y_  » 0 plane  at  the  point  where  it  is  intersected  by  the 
line  connecting  the  leading  edges  of  the  root  chords  of  the  left  and 
right  wing  halves.  The  coordinate  system  (x,y,z)  with  its  origin  at 
corner  1 of  the  cross-hatched  panel  is  associated  with  a semi-infinite 
triangle  which  is  part  of  the  solution  for  the  panel  to  be  described 
later. 

A distribution  of  constant  u-velocity  type  panels  to  model  loading 
and  a distribution  of  constant  source  type  panels  to  model  thickness  will 
be  laid  out  on  the  wing  and  pylon.  The  distributions  of  the  two  types 
of  panels  need  not  be  the  same  in  the  chordwise  direction,  but  they  do 
coincide  in  the  spanwise  direction.  The  sides  of  both  types  of  panels 
are  taken  parallel  to  the  centerline  of  the  fuselage.  Their  leading  and 
trailing  edges  are  swept  along  constant  percent  chordlines  of  the  wing 
or  pylon.  A leading  or  trailing  edge  is  called  subsonic  or  supersonic 
depending  on  whether  or  not  the  component  of  the  tree-stream  velocity 
perpendicular  to  the  edge  in  question  is  subsonic  or  supersonic, 
respectively.  Figure  1 shows  a simplified  layout  of  2 chordwise  by 
b spanwise  panels  on  the  left  wing  half.  A distribution  of  2 chordwise 
ay  2 spanwise  panels  covers  the  pylon.  Generally,  more  panels  are 
required  to  obtain  accurate  results  as  described  in  the  section  concerned 
with  comparisons. 

Constant  u-velocity  panels  only  are  laid  out  over  part  of  the 
fuselage  and  designated  body  interference  panels.  The  leading  and 
trailing  edges  of  these  panels  are  unswept.  Only  the  left  half  of  the 
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aircraft  configuration  needs  to  be  covered  because  of  symmetry  about 
the  y « 0 plane  provided  that  the  effects  of  the  other  half  are 
accounted  for  as  discussed  later.  The  length  of  fuselage  covered  by 
the  body  interference  panels  is  determined  by  the  range  of  influence  of 
the  wing,  pylon,  and  store.  This  range  is  taken  as  the  length  along  the 
body  measured  from  the  leading  edge  of  the  wing-fuselage  junction,  or 
root  chord,  to  the  trailing  edge  for  a supersonic  wing  trailing  edge. 

For  the  case  of  a subsonic  trailing  edge,  the  length  should  extend  back 
further.  In  this  case,  the  interference  shell  should  be  extended  back 
as  far  as  the  axial  location  of  the  trailing  edge  of  the  wing  tip.  As 
such,  the  length  of  the  body  interference  shell  is  sufficient  to  account 
for  wing-pylon  to  fuselage  interference.  If  the  store  to  be  separated 
lies  close  to  the  fuselage,  the  fuselage  interference  shell  must  also 
be  long  enough  to  cancel  store  effects  on  the  fuselage.  In  the  simplified 
layout  of  figure  1,  the  body  interference  shell  consists  of  5 circum- 
ferential rings  with  4 pan*»’d  per  half  ring. 

Each  constant  u-velocity  panel  contains  a control  point  placed  at 
the  95-percent  location  on  the  chord  containing  the  panel  centroid  as 
shown  in  figure  1.  The  flow  tangency  condition  is  applied  at  these 
points,  figure  2 shows  the  geometric  characteristics  of  a panel  in 
detail.  Before  specifying  the  boundary  condition,  the  method  of  obtaining 
the  constant  u-velocity  and  constant  source  solutions  for  one  panel  will 
be  described. 

3.3.2  Constant  u-velocity  and  constant  source  panel  solutions 

An  explanation  of  the  basic  solution  for  the  two  types  of  trapezoidal 
panels  is  now  given  for  completeness.  The  treatment  is  a condensed  version 
of  the  one  given  in  reference  10.  However,  all  the  equations  necessary 
for  the  solutions  are  contained  in  this  report. 

In  reference  9,  Woodward,  et  al . , derived  ganeral  expressions  for 
the  solution  to  the  wave  equation  for  a semi-infinite  triangular  shape 
with  a sweptback  leading  edge  subject  to  a jump  condition  in  axial  (or  u) 
velocity.  The  reference  also  contains  the  solution  for  the  semi-infinite 
triangle  with  constant  source  strength.  As  written  by  Woodward,  these 
results  have  a removable  singularity  for  zero  leading-edge  sweep.  The 
method  of  deriving  these  expressions  is  described  in  Appendix  II,  and 
the  forms  of  the  expressions  used  in  the  computer  programs  are  given. 

The  layout  of  the  semi-infinite  triangle  with  either  a sweptback  or 


sweptforward  leading  edge  and  its  relationship  to  the  free-stream  Mach 
cone  are  shown  in  figure  3.  The  sweptback  case  shown  corresponds  to  a 
subsonic  leading  edge.  The  coordinate  system  associated  with  the  tri- 
angle is  also  shown.  Note  that  the  x and  z directions  are  opposite  to 
the  xw  and  zw  directions  shown  in  figure  1.  The  triangle  lies  in  the 
z m o or  planform  plane.  It  should  also  be  noted  that  the  side  edge 
lies  along  the  x axis.  The  available  solutions  in  Appendix  II  pertain 
to  the  3weptback  triangular  shape. 

The  solution  for  a constant  u-velocity  or  constant  source  panel, 
such  as  the  cross-hatched  one  shown  in  figure  1,  is  obtained  through  a 
superposition  scheme  using  four  semi-infinite  triangular  shapes  with 
their  apexes  at  each  of  the  corners.  The  superposition  scheme  applied 
to  the  cross-hatched  panel  in  figure  1 with  sweptback  leading  and  trailing 
edges  on  the  left  wing  is  shown  in  figure  4(a) . The  coordinate  system 
shown  in  figure  4(a)  is  consistent  with  the  coordinate  system  associated 
with  the  serai-infinite  triangle  of 'figure  3,  but  the  orientations  of  the 
semi-infinite  triangles  placed  at  the  corner  points  correspond  to  swept- 
forward triangles.  Thus,  a transformation  is  used  which  reverses  the  y 
orientation  of  the  triangles.  The  available  solutions  in  Appendix  II 
are  then  applied  and  the  sign  on  the  resulting  sidewash  is  changed.  If 
the  panel  leading  or  trailing  edge  is  sweptforward  (as  on  the  right  wing) , 
the  scheme  such  as  that  depicted  in  figure  4(b)  for  that  edge  must  be 
employed  in  the  superposition  method  for  a trapezoidal  panel  on  the  left 
wing.  In  this  case,  the  required  solution  for  each  semi-infinite  triangle 
corresponds  directly  to  the  triangle  with  sweptback  leading  edge  shown  in 
figure  3 and  is  given  in  Appendix  II.  Trapezoidal  panels  may  have  swept- 
back, sweptforward,  or  mixed  leading  and  trailing  edges.  Therefore,  the 
superposition  methods  used  in  the  computer  program  handle  these  conditions 

The  superposition  principle  can  also  be  used  to  obtain  the  pertur- 
bation velocities  due  to  a constant  u-velocity  or  constant  source  panel 
at  any  point  in  the  field.  Woodward,  et  al.,  in  reference  9 give 
expressions  for  the  velocity  components  due  to  a semi-infinite  triangle 
with  sweptback  leading  edge  which  can  be  used  in  the  superposition 
process.  The  alternate  forms  of  these  expressions  used  in  the  computer 
program  are  given  in  Appendix  II. 

As  described  in  Appendix  II,  the  source  panel  strengths  are  directly 
related  to  the  streamwise  slopes  of  the  wing  and  pylon  thickness  distri- 
butions and  therefore  must  be  specified.  For  blunt  leading  edges,  a 
special  procedure  (described  later)  is  used. 


With  the  induced  velocity  expressions  for  the  two  types  of  panels 
known,  it  is  possible  to  formulate  the  influence  of  a constant  u-velocity 
panel  of  unknown  strength  and  the  influence  of  a constant  source  panel 
with  known  strength  at  any  control  point  accounting  for  the  region  of 
influence  inherent  in  supersonic  flight.  By  considering  interactions 
between  all  constant  u-velocity  panels  laid  out  on  the  wing .pylon, and 
body  interference  shell,  there  results  a set  of  simultaneous  equations 
from  which  the  unknown  constant  u-velocity  panel  strengths  can  be 
determined.  This  formulation  will  now  be  discussed. 

3.3.3  Flow  tangency  boundary  condition 

The  flow  tangency  boundary  condition  states  that  there  is  no  flow 
through  the  constant  u-velocity  panels  on  the  wing, pylon,  and  interference 
shell  at  each  control  point.  A few  control  points  are  shown  in  figure  1. 

Velocities  normal  to  the  wing  include  a component  of  the  free-stream 
and  perturbation  velocities,  ww,  induced  by  the  constant  u-velocity 
panels  distributed  over  the  wing,  pylon,  and  interference  shell.  The 
line  sources,  sinks  and  line  doublets  used  to  model  the  fuselage  induce 
interference  velocities  at  the  control  points  distributed  over  the  wing. 
Components  of  these  velocities  normal  to  the  wing  are  also  included  in 
the  wing  boundary  condition.  In  this  way,  fuselage -to -wing  interference 
is  accounted  for.  Furthermore,  the  influence  of  the  store  body  model  is 
included  in  the  wing  boundary  condition.  Lastly,  the  effects  of  the 
source  panels  laid  out  on  the  pylon  are  added  to  the  wing  boundary  con- 
dition to  account  for  pylon  thickness  interference  on  the  wing.  The 
velocities  induced  by  the  fuselage  and  store  sources,  sinks,  and  doublets 
and  wing  thickness  are  designated  externally  induced  interference  down- 
wash  w^..  Let  v be  the  control  point  index  for  the  constant  u-velocity 
panels  on  the  wing  and  let  n he  the  summation  index  for  all  the 
influencing  constant  u-velocity  panels.  The  boundary  condition  applied 
at  control  points  on  the  wing  is  then  given  in  terms  of  the  wing  coordi- 
nate system  (x^,y^,z  ) shown  in  figure  1 as  follows 


NPTOT 
n-l 

v - 1,2, . . .N PAULS 
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where  NPTOT  is  the  sum  of  the  number  of  constant  u-velocity  panels  on 
the  wing,  pylon,  and  fuselage  interference  shell.  The  number  of  such 
panels  on  the  wing  is  NPANLS.  Angle  af  is  the  nondimens ional  upwash 
due  to  the  free  stream  at  the  wing  chordal  plane  and  angle  is  the 

local  angle  of  attack  due  to  wing  camber  and/or  twist  as  shown  in  the 
following  sketch. 


The  boundary  condition  on  the  pylon  states  that  the  net  velocity 
through  the  constant  u-velocity  panels  at  their  control  points  equals 
zero.  Velocities  normal  to  the  pylon  include  perturbation  velocities 
v induced  by  the  constant  u-velocity  panels  distributed  over  the  wing, 
pylon,  and  interference  shell.  Components  of  perturbation  velocities 
induced  by  the  line  sources,  sinks  and  line  doublets  representing  the 
fuselage  are  included.  Thus,  fuselage-to-pylon  interference  is  accounted 
for.  Effects  of  the  store  line  sources,  sink3  and  line  doublets  are 
also  included.  The  pylon  boundary  condition  also  contains  contributions 
from  the  source  panels  on  the  wing  accounting  for  wing  thickness  inter- 
ference on  the  pylon.  The  last  three  contributions  to  sidewash  added 
together  are  called  externally  induced  sidewash  vWj_-  Taking  v as 
| the  index  of  the  control  point  on  the  pylon  and  n as  the  index  of  the 

i influencing  panels,  the  boundary  condition  is  expressed  as 


NPTOT  v. 
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n=l 


W 


v.n 


v. 
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Wi.v 

v_ 


(10) 


v ■=  NPANLS+1 , NPANLS-t-2 NPANLS+MP 

This  condition  is  also  expressed  in  the  wing  coordinate  system  of  fig- 
ure 1.  The  number  of  constant  u-velocity  panels  distributed  on  the  pylon 
is  MP.  The  right-hand  side  of  equation  (9)  represents  the  free-stream 
component  and  externally  induced  perturbation  downwash  and  the  right-hand 
side  of  equation  (10)  represents  externally  induced  perturbation  sidewash. 
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Similarly,  letting  v^v  n be  the  perturbation  velocity  in  the  direction 
normal  to  the  body  interference  panel,  the  flow  tangency  condition 

states  that  the  net  velocity  normal  to  the  body  interference  panel  is 
zero. 
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V - NPANLS+MP+l,NPANIiS+MP+2,  . . .NPTOT 

(i  **  v-NPANLS-MP 


The  right-hand  side  of  the  above  equation  represents  the  externally 
induced  perturbation  velocity  normal  to  the  body  interference  panel 
under  consideration.  It  contains  interference  velocities  induced  by 
source  panels  laid  out  on  the  wing  and  the  pylon  to  model  thickness. 
Store  line  sources,  sinks  and  line  doublet  effects  are  also  included. 
In  this  way,  interference  on  the  fuselage  due  to  wing-pylon  thickness 
and  store  body  are  accounted  for.  Angle  0 is  called  the  body  panel 
orientation  angle  measured  relative  to  a plane  parallel  to  the  (xB,yB) 
or  (x^y^  plane  as  indicated  in  the  following  sketch  looking  in  the 
positive  Xg  direction.  The  angle  0 shown  in  the  sketch  is  a 
negative  angle.  The  summation  on  the  left-hand  side  of  equation  (9) 


represents  the  downwash  perturbation  velocity  at  the  control  points  on 
the  wing  induced  by  the  constant  u-velocity  panels  on  the  wing,  pylon, 
and  interference  shell.  The  summation  on  the  left-hand  side  of  equa- 
tion (10)  represents  the  sidewash  perturbation  velocity  at  the  control 
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points  on  the  pylon  induced  by  all  the  constant  u-velocity  panels. 

Finally,  the  summation  on  the  left-hand  side  of  equation  (11)  represents 
the  perturbation  velocity  normal  to  the  interference  panel  under  consid- 
eration induced  by  all  constant  u-velocity  panels.  So  far,  the  strengths 
of  these  panels  are  still  unknown. 

in  the  boundary  conditions  formulated  above,  advantage  is  taken  of 
the  fact  that  the  yB  • 0 or  y^  « 0 plane  is  a plane  of  symmetry.  Panels 
are  laid  out  over  the  left  wing  and  pylon  and  only  the  left  half  of  the 
fuselage  is  covered  with  interference  panels.  Perturbation  velocities 
on  both  sides  of  equations  (9) , (10) , and  (11)  are  induced  not  only  by 
the  panels  to  the  left  of  the  symmetry  plane  but  also  should  contain 
contributions  from  the  right  half.  The  required  procedure  is  discussed 
below  in  connection  with  the  aerodynamic  influence  coefficients. 

Next,  it  will  be  shown  how  the  perturbation  velocities  are  expressed 
in  terms  of  aerodynamic  influence  coefficients.  These  coefficients 
relate  the  perturbation  velocity  components  induced  at  some  point  by  a 
constant  u-velocity  panel  to  its  strength  u+/VM  and  the  coordinates  of 
the  point  relative  to  the  panel  corners.  A similar  procedure  relates  the 
perturbation  velocities  at  a field  point  to  the  known  strength  of  a source 
panel . 

3.3.4  Aerodynamic  influence  coefficients 

Appendix  II  contains  the  results  of  the  theory  for  determining  the 
pert  irbation  velocities  induced  by  a semi- inf inite  triangle  subject  to 
a constant  jump  condition  in  axial  velocity.  The  appendix  also  specifies 
the  perturbation  velocities  induced  by  a semi-infinite  triangle  with  a 
constant  strength  source  distribution.  The  results  are  given  below  in 
symbolic  form.  In  terms  of  the  coordinate  system  associated  with  the 
triangle  shown  in  figure  3,  the  perturbation  velocities  due  to  a triangle 
with  sweptback  leading  edge  and  with  constant  jump  in  u-velocity  are 
expressed  as  follows. 


V!o  “ TT  (v^)  FU(X'y’Z’0’’W 
\T  ” •ii:(iV^)Fw(x’y‘z»0*^LE) 
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The  quantity  u+/VK  is  the  constant  strength  of  the  singularity  distri- 
buted over  the  triangle.  In  fact,  it  is  proportional  to  the  axial 
velocity  u+  in  the  plane  of  the  triangle  (2  ««  0) . The  computer  program 
has  been  arranged  to  predict  axial  velocity  u+  in  the  zero  plane  and 
requires  a very  small  negative  value  for  2 to  produce  axial  velocity 
u shown  in  figure  3.  The  complete  expressions  for  the  perturbation 

velocities  are  given  in  Appendix  II.  Functions  F.  ,F  ,F  are  called 

u v w 

influence  functions  and  depend  on  coordinates  x,y,z  of  the  point  at 
which  the  velocities  are  computed  relative  to  the  apex  of  the  semi- 
infinite triangle  shown  in  figure  3.  The  influence  functions  can  be 
determined  from  Appendix  II,  equations  (II-4)  and  (11-12),  and  depend 
on  the  leading-edge  sweep  il/T  „ and  the  factor  3 given  in  terms  of 
the  free-stream  Mach  number  by 


V M=  - 1 


(13) 
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Perturbation  velocities  induced  by  the  semi-infinite  triangle  with 
constant  source  strength  are  related  to  a surface  tangent,  tan  6,  to  be 
discussed  later,  and  the  coordinates  of  the  field  point  relative  to  the 
triangle  apex.  Employing  the  subscript  t to  indicate  association  with 
thickness,  the  perturbation  velocities  are  expressed  as 
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Functions  Fu^.,FVt,FWt  are  the  influence  functions  associated  with  the 
thickness  and  depend  on  coordinates  x,y,z  of  the  field  point  relative 
to  the  apex  of  the  semi-infinite  triangle  such  as  the  one  with  sweptback 
leading  edge  shown  in  figure  3.  These  functions  can  be  determined  from 
Appendix  II,  equations  (n-15)  and  (II-l6),and  depend  on  the  leading-edge 
sweep  ^LE  and  the  factor  p given  by  equation  (13) . 
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In  order  to  formulate  the  perturbation  velocities  induced  by  a 
trapezoidal  constant  u-velocity  panel,  the  superposition  principle 
described  earlier  and  depicted  in  figure  4 must  be  employed.  For  example, 
for  a panel  with  sweptback  leading  and  trailing  edges  on  the  left  wing, 
the  upwash  at  the  control  point  induced  by  the  n^  constant 

u-velocity  panel  is  expressed  in  accordance  with  figure  4(a)  as 


- Fwn(xi , v ’ ”^i , v ,zi  ,v 

Fw  ^X4,v’  '^4JV*Z41V*^*^0KE^  + FW  ^X3,V*  ^3 , V * Z 3 , V ’ ^ * ^TE^  J 


The  subscripted  coordinates  represent  the  coordinates  of  the  control 
point  relative  to  one  of  the  four  panel  corner  points.  The  numbering 
sequence  of  the  corner  points  is  kept  consistent  in  its  relation  to  the 
directions  of  the  coordinate  systems  shown  in  figure  4.  The  leading-edge 
and  trailing-edge  sweep  angles  of  the  panel  are  designated  and 

respectively.  They  are  indicated  in  figure  2.  Here  function  Fw  is 
defined  in  equation  (12) . Similar  expressions  can  be  written  for  u/V^ 
and  -v/VM  using  Fu  and  Fv,  respectively.  Likewise,  the  perturbation 
velocities  induced  by  a trapezoidal  constant  source  panel  are  obtained 
by  means  of  a superposition  scheme  using  the  solutions  for  four  semi- 
infinite triangles  with  their  apexes  at  each  of  the  panel  corners.  The 
source  strength  for  one  panel  is  given  by  the  surface  slope  of  the  wing 
or  pylon  thickness  distribution  at  the  panel  centroid.  It  is  expressed 
as  the  tangent  of  the  angle  9 associated  with  the  thickness  envelope 
as  shown  in  the  following  sketch.  Therefore,  the  upwash  due  to  thickness 


Thickness  envelope 


Mean  camber  surface 
made  to  be  flat 


Tangent  to  thickness  envelope 
at  axial  location  of  the 
centroid 


Source  panel  centroid 


at  the  vth  control  point  induced  by  the  constant  source  panel  with 
index  nt  on  the  left  wing  is  expressed  in  accordance  with  figure  4(a) 
as 


Fwt,nt('xi)v,’’yiJv,:zi,v,i3,^tLE^  " Fwt>nt^X4,v,-ir4,v,:i4,v*P,^tTB^ 

V„t<*3,v>-r0,v>23,v'e.%E']  <16) 


Here  pertains  to  the  leading-edge  or  trailing-edge  sweep  of  the 

constant  source  panel  accounting  for  thickness.  Note  that  if  either  the 
leading  or  trailing  panel  edge  is  sweptf orward , the  scheme  shown  in 
figure  4(b)  must  be  used  for  the  edge  in  question.  The  bracketed  terms 
in  equations  (15)  and  (16)  are  designated  as  the  upwash  aerodynamic 
influence  coefficients  associated  with  a constant  u-velocity  panel  and  a 
constant  source  panel,  respectively. 

3.3.5  Symmetry  considerations 

In  order  to  account  for  the  effects  of  the  constant  u-velocity  and 
source  panels  on  the  right  half  of  the  aircraft  configuration,  the  per- 
turbation velocity  components  induced  by  a panel  on  the  left  are  augmented 
by  a contribution  from  its  image  panel  on  the  right.  However,  the 
following  simplification  is  used.  The  effect  of  the  image  panel  at  a 
given  point  can  be  obtained  using  the  panel  on  the  left  with  a change 
in  sign  of  the  yw  or  yg  coordinate  of  the  point.  The  sidewash  calcu- 
lated this  way  at  (xw, -yw, zw)  is  then  transformed  back  by  a change 

in  sign.  The  procedure  is  shown  in  the  following  sketch  for  a panel  on 
the  left  wing  half  of  a wing-alone  configuration.  The  coordinate  system 
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(x,y,z)  associated  with  the  semi-infinite  triangle  with  its  apex  at 
corner  1 will  be  used.  It.  is  the  same  as  the  coordinate  system  of  fig- 
ure 3.  The  coordinates  of  point  P relative  to  corner  1 ,yw  ,zw  ) 
are  then  given  by  l i i 


X - -(xw  - x^) 

y - % - yVi  > 


(17) 


With  these  local  coordinates,  the  direct  influence  of  corner  1 at  point 
P can  be  calculated.  The  influences  of  the  other  corners  of  the  panel 
are  added  in  accordance  with  the  superposition  scheme  as  indicated  for 
the  w component  in  equations  (15)  and  (16) . The  coordinates  of  the 
point  (xw,-yw,zw)  relative  to  corner  1 are 

- -(xw  - *W  0 
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The  influence  of  corner  1 ia  now  calculated  at  point  P^.  The  influences 
of  the  other  corners  of  the  panel  are  added  in  accordance  with  the  super- 
position scheme.  The  resulting  velocity  components  u,v,w  represent 
the  influence  of  the  panel  on  the  left  wing-half  on  the  point  with 

a change  in  the  sign  of  the  sidewash  v,  they  also  represent  the  influence 
of  the  image  panel  on  the  point  P.  Thus,  at  a given  point  P there 
will  be  a direct  influence  induced  by  a constant  u-velocity  or  constant 
source  panel  and  an  indirect  effect  induced  by  the  image  of  the  panel 
under  consideration. 

This  method  of  accounting  for  symmetry  is  also  applied  to  constant 
u-velocity  and  constant  source  panels  on  a wing-mounted  pylon  and  the 
constant  u-velocity  panels  on  the  interference  shell.  The  influence 
functions  required  in  expressions  (15)  and  (16)  and  defined  by  equa- 
tions (12)  and  (14)  for  these  panels  are  obtained  from  Appendix  II,  equa- 
tions (II-4) , (11-12)  and  equations  (11-15) , (11-16) , using  suitable 
transformations.  For  the  pylon,  a 90°  counterclockwise  rotation  about 
the  x-axis  (viewing  upstream)  places  the  semi-infinite  triangles  in  the 
vertical  plane  and  the  influence  functions  can  be  determined.  The  super- 
position scheme  gives  the  solution  for  a trapezoidal  panel  on  the  pylon, 
affects  of  the  image  wing-mounted  pylon  panel  are  obtained  in  the  manner 
described  above  for  a panel  on  the  wing.  The  same  procedure  is  applied 
to  the  pylon  source  x^anels.  The  influence  functions  associated  with  a 
body  interference  panel  are  also  obtained  after  performing  a rotation 
about  the  x-axis.  The  angle  of  rotation  is  different  ft r each  panel  and 
is  related  to  the  body  orientation  angle  0^  mentioned  .arlier  in  connec- 
tion with  the  boundary  conditions.  Image  panel  effects  are  accounted  for 
using  the  same  procedure  as  used  for  a panel  on  the  wing. 

The  terms  on  both  sides  of  equations  (9)  , (10)  , and  (11)  can  nov; 
be  specified.  Equation  (15)  for  w/VM  and  a similar  expression  for 
-v/V  n.:e  used  on  the  left-hand  sides.  A combination  of  the  twe  compo- 
nents is  required  in  equation  (11) . On  the  right-hand  side,  the  contri- 
butions from  the  fuselage  line  sources,  3inks  and  line  doublets  are  given 
by  equations  (1-30)  in  Appendix  I.  store  effects  are  given  by  the  same 
equations  but  require  a translation  of  coordinate  axes.  Wing  and  pylon 
source  panel  effects  are  given  by  equation  (16)  for  w^/V^  with  similar 
expressions  for  u^/V^  and  -vt/V(n.  Because  of  the  difference  in  coordi- 
nate systems  shown  in  figures  1 and  3,  the  following  transformation  is 
required  to  obtain  the  perturbation  velocities  in  the  wing  coordinate  system. 
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Perturbation  velocities  associated  with  the  constant  u-velocity  panels 
must  also  be  transformed  as  follows  into  the  wing  coordinate  system. 
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3.3.6  Solution  of  constant  u-velocity  panel  strengths 

After  recasting  the  perturbation  velocities  on  both  sides  of  equa- 
tions (9) , (10) , and  (11)  in  terms  of  the  influence  functions  specified 
by  equations  (12)  and  (14)  accounting  for  symmetry  and  the  transformations 
in  equations  (19)  and  (20),  there  results  a set  of  simultaneous  equations 
in  which  the  unknowns  are  the  NPTOT  values  of  panel  strengths  u+  A„. 

The  values  can  be  obtained  through  a matrix  solution  for  given  angle  of 

attack,  a,,  and  Mach  number,  M'  . 

1 f '00 

3.4  Flow  Field  Calculation  Including  the  Effects  of  the  Shock 
Associated  with  the  Wing  Leading  Edge 

The  method  used  to  calculate  the  flow  field  below  the  wing  plane  of 
the  parent  aircraft  is  discussed  in  this  section.  The  flow  models  repre- 
senting the  wing,  pylon,  and  fuselage  in  supersonic  flow  are  based  on 
linear,  potential  flow  theory  and  are  described  in  the  preceding  sections. 
In  order  to  predict  the  flow  field  in  the  vicinity  of  the  airframe  with 
some  degree  of  accuracy,  the  linear  flow  theory  must  be  augmented  by 
nonlinear  corrections.  In  linear  theory,  all  disturbances  travel  and 
attenuate  along  Mach  waves  in  the  two-dimensional  case  or  Mach  cones  in 
the  three-dimensional  case.  For  infinitely  small  disturbances,  the 
description  of  the  flow  field  and  weak  shock  waves  by  linear  theory  is 
sufficiently  accurate.  However,  when  the  disturbances  are  large,  they 
will  travel  along  curved  characteristics  in  accordance  with  nonlinear 
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theory.  Detached  and/or  attached  shocks  associated  with  the  fuselage 
nose  and  wing  leading  edge  will  not  lie  along  Mach  lines  or  Mach  cones. 

In  this  work,  a simplified  method  is  used  to  position  the  wing 
leading-edge  shock  only.  The  technique  is  based  on  the  assumption  that 
the  shock  pattern  from  the  leading  edge  depends  on  wing  thickness  only. 

For  nonzero  angle  of  attack,  the  shock  pattern  is  assumed  to  rotate  with 
the  wing.  This  assumption  is  based  on  experimental  evidence  discussed 
later  in  the  comparisons. 

For  a wing  with  either  a blunt  or  sharp  leading  edge,  the  flow  field 
due  to  wing  thickness  only  is  computed  along  a traverse  parallel  to  the 
wing  chord  at  the  required  ywizw  location.  The  flow  properties  at  the 
location  whore  wing  effects  due  to  thickness  are  first  felt  are  used  to 
determine  a local  Mach  number.  The  flow  properties  due  to  wing  thickness 
along  the  traverse  are  then  recalculated  with  the  local  Mach  number.  The 
location  where  wing  thickness  effects  are  now  felt  first  is  taken  as  the 
location  of  the  shock.  This  procedure  does  not  depend  on  the  wing  leading 
edges  being  subsonic  or  supersonic.  The  procedure  used  to  determine  the 
flow  field  induced  by  the  entire  airframe  will  now  be  summarized.  A 
detailed  description  is  given  in  the  next  section.  The  effects  of  wing 
thickness  and  lift  are  computed  at  the  field  point  in  question  on  the 
basis  of  the  local  Mach  number.  The  contribution  from  the  interference 
shell  is  also  calculated  using  the  local  Mach  number.  Velocities  induced 
by  the  fuselage  sources,  sinks  and  doublets  are  still  calculated  using 
the  free-stream  Mach  number.  If  a pylon  is  present,  its  effects  are  also 
based  on  the  free-stream  Mach  number. 

3.4.1  Flow  angle  and  local  Mach  number  calculation 

Consider  the  basic  airframe  configuration  comprised  of  a wing  with 
thickness  attached  to  a fuselage.  The  first  step  is  to  solve  for  the 
singularity  strengths  representing  this  configuration  in  accordance  with 
the  methods  based  on  linear  theory  described  earlier  for  the  flight 
conditions  at  hand.  The  result  is  a complete  flow  model  made  up  of  line 
sources,  sinks  and  line  doublets  to  model  the  fuselage  and  constant  u- 
velocity  panels  to  model  the  lifting  wing  and  wing-fuselage  interference. 
There  will  also  be  source  panels  to  account  for  wing  thickness.  The 
strengths  of  all  the  singularities  are  known  at  this  point  and  the 
complete  flow  model  is  based  on  linear,  potential  flow  theory.  Special 
care  must  be  exercised  when  assigning  values  to  the  strengths  of  the 
wing  source  panels  as  described  next. 
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The  strength  of  each  source  panel  is  directly  related  to  the  slope, 
tan  Q,  of  the  wing  thickness  envelope  evaluated  at  the  panel  centroid. 

For  wings  with  supersonic,  sharp  leading  edges  for  which  the  shock  remains 
attached,  the  magnitudes  of  the  slopes  at  the  leading  edge  are  bounded. 
Figure  4 in  reference  11  contains  a curve  relating  the  minimum  Mach  number 
to  wedge  angle  for  the  shock  to  be  attached.  Since  the  leading  edge  of 
the  wing  is  supersonic,  the  pertinent  Mach  number  and  wedge  angle  are 
defined  in  the  direction  perpendicular  to  the  wing  leading  edge.  They 
are  related  to  the  streamwise  properties  as  follows. 
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Here  is  the  sweep  angle  of  the  wing  leading  edge.  If  for  the 

given  free-stream  Mach  number  and  wing  thickness  slope  at  the  leading 
edge  the  shock  is  detached,  the  procedure  described  below  is  used.  This 
procedure  is  also  used  for  all  wings  with  supersonic,  blunt  leading  edges 
and  involves  certain  concepts  from  normal  shock  analysis  summarized  in 
Appendix  111.  The  maximum  slope  input  in  the  computer  program  is  related 
to  the  shock  detachment  wedge  angle,  fidet.  This  angle,  ddut,  is  the 
semi-vertex  angle  of  a wedge  for  which  the  shock  first  detaches  for  a 
given  Mach  number  normal  to  the  leading  edge.  Appendix  III  contains  a 
graph  baaed  on  chart  2 of  reference  11  from  which  6det  can  be  deduced. 
The  maximum  slope  tan  8max  in  the  streamwise  direction  is  then  given  by 


tan  Gmax  " (tan  6det)coa  ^WLE  <22 

If  for  a given  wing  the  thickness  slope  of  one  of  the  forward  source 
panels  is  larger  than  tan  0m,  . the  slope  is  set  equal  to  tan  0ma„ 
in  the  input  data. 

The  case  of  the  subsonic  leading  edge  will  now  be  discussed.  In 
this  case,  the  flow  in  the  direction  normal  to  the  wing  leading  edge  is 
subsonic.  Considering  only  linear  theory,  there  is  a limiting  Mach  cone 
(based  on  free-stream  Mach  number)  with  its  origin  at  the  leading  edge 
of  the  wing-fuselage  junction  as  shown  in  the  following  sketch.  On  the 
basis  of  a local  Mach  number  discussed  later,  the  linear  theory  is 
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j corrected.  The  corrected  theory  predicts  nonzero  influences  ahead  of  the 

limiting  Mach  wave.  Thus,  in  a nonlinear  sense,  a shock  shape  is  gene- 
rated which  lies  in  front  of  the  limiting  Mach  cone.  In  fact,  this  shock 
is  the  envelope  of  many  nonlinear  Mach  "cones." 

If  the  subsonic  leading  edge  of  a wing  is  sharp,  the  maximum  stream- 
wise  slope,  tan  0(nax,  is  the  value  at  the  leading  edge.  If  the  leading 
edge  is  blunt,  the  thickness  slopes  input  into  the  computer  program  must 
be  bounded.  Since  the  Mach  number  normal  to  the  leading  edge  is  sub- 
sonic, the  method  given  above  for  the  supersonic  lending  edge  is  no 
longer  valid.  However,  assuming  that  the  shock  or  ginates  at  or  near 
the  leading  edge  of  the  wing-fuselage  junction,  an  approach  analogous 
to  the  one  used  for  the  supersonic  leading-edge  case  is  employed.  The 
supersonic  flow  runs  parallel  to  the  fuselage  as  it  approaches  the 
leading  edge  of  the  wing-fuselage  junction.  The  shock  detachment  angle, 
^det’  determined  for  the  streamwise  Mach  number,  and  the  maximum 
slope,  tan  0max,  is  then  taken  equal  to  tan  $det  directly.  Away  from 
the  junction,  this  upper  bound  is  retained  for  all  spanwise  locations. 

If  for  a given  blunt  wing  with  a subsonic  leading  edge  the  streamwise 
thickness  slope  of  one  of  the  forward  source  panels  is  larger  than 

tan  0 , the  slope  is  set  equal  to  tan  0 

max  ^ max 

; For  a given  field  point  (xw,yw,zw)  under  the  wing  at  which  the  flow 

field  is  to  be  calculated,  an  axial  traverse  is  laid  out  parallel  to  the 
fuselage  centerline  (and  therefore  parallel  to  the  wing  plane) . in  the 

I 

crossflow  plane,  the  lateral  coordinates  of  the  traverse  are  (yT,,zr.) 

JW’  W' 

associated  with  the  field  point  under  consideration.  For  the  purpose  of 
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calculating  the  local  Mach  number,  the  length  of  the  traverse  is  chosen 
to  be  from  the  axial  location  of  the  leading  edge  of  the  wing-fuselage 
junction  (or  wing  root  chord)  to  the  axial  location  of  the  trailing  edge 
of  the  local  wing  chord.  The  local  wing  chord  is  the  chord  which  has 
the  lateral  coordinates  (yw,0)  in  the  crossflow  plane.  For  a sharp  or 
blunt  wing  with  either  a subsonic  or  supersonic  leading  edge,  the  flow 
field  due  to  wing  thickness  only  is  computed  at  points  along  the  traverse 
from  its  starting  point  to  its  end  by  means  of  the  linear  theory  method 
described  in  an  earlier  section.  For  example,  the  upwash  due  to  one 
source  panel  is  given  by  equation  (16) . All  three  velocity  components 
due  to  wing  thickness  are  then  obtained  by  summing  the  effects  of  all 
the  source  panels  laid  out  over  the  left  wing  and  their  images  on  the 
right  wing. 

Once  the  perturbation  velocities  due  to  wing  thickness  only  have 
been  calculated  along  the  traverse,  the  point  on  the  traverse  at  which 
the  effects  are  first  felt  is  isolated.  This  point  corresponds  to  the 
intersection  of  the  traverse  with  a Mach  cone  based  on  frae-stream  Mach 
number  and  with  its  origin  at  one  corner  of  a source  panel.  This  corner 
lies  on  the  leading  edge  and  inboard  of  the  traverse  for  a sweptback  wing 
leading  edge.  The  situation  is  depicted  below  in  side  view. 
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The  change  in  flow  angle,  Av,  is  determined  at  the  point  on  the 
traverse  where  wing  thickness  influence  is  first  felt.  The  quantity 
Av  can  be  calculated  from  the  velocity  components,  u^/V^,  and 

vv„- 
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tan  Av 


(23) 


Velocity  component  u/VM  is  the  nondimens ional  total  velocity  in  the 
free-stream  direction  and  and  v^/V^  are  the  lateral  perturbation 

velocities  in  the  wing  coordinate  system  shown  in  figure  1.  In  this 
system,  the  nondimensional  total  axial  velocity  is  then  given  by 
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The  change  in  flow  angle,  Av,  corresponds  to  a compression  of  the  free- 
streara  flow.  Assuming  that  this  compression  is  isentropic,  the  Prandtl- 
Meyer  angle  at  the  point  where  the  flow  is  first  disturbed  is  obtained 
from  Prandtl-Meyer  flow  theory. 

v2  - v - Av  (25) 

Prandtl-Meyer  angle  v2  is  related  to  the  free-stream  Mach  number 

M (”  M ) . It  is  given  by  equation  (171c)  in  reference  11  and  repeated 
00  r 1 

here  for  convenience. 

Vi  - 2.4495  tan"1 [o. 40825  (M*  - 1) 1/2  ] - tan"1 (M2  - 1) l/s  (26) 


In  accordance  with  reference  12,  the  Mach  number  associated  with  Prandtl- 
Meyer  angle  vg  can  be  determined  most  conveniently  as  follows  in 
accordance  with  the  Prandtl-Meyer  relationship. 
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1 + 1.3604  p + 0.0962  p2  - 0.5127  p3 
1 - 0.6722  p - 0.3278  p2 
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The  author  of  reference  12  claims  0.05-percent  accuracy  over  the  entire 
range  of  Mach  number  (1  < Mg  < ») . The  component  of  the  Mach  number, 
M2,  in  the  axial  or  free-stream  direction  is  then  given  by 


26 


(28) 


M2  - M£cos  Av  “ M . 

streamwise  1 

From  here  on,  this  streamwise  Mach  number  will  be  designated  and  referred 
to  as  the  local  Mach  number,  M^. 

3.4.2  Shock  shape  and  flow  field  calculated  on  the  basis  of 
local  Mach  number 

Under  the  assumption  that  the  shape  of  the  shock  associated  with  the 
wing  leading  edge  depends  on  wing  thickness  only,  the  scheme  to  be 
described  next  allows  for  the  determination  of  the  shock  location.  The 
assumption  will  be  justified  by  evidence  offered  by  comparison  with 
experimental  data  which  will  be  discussed  later. 

The  summed  effects  of  the  source  panels  on  the  wing  are  computed 
again  at  points  along  the  same  axial  traverse  discussed  previously. 
However,  the  Mach  number  used  in  this  calculation  is  calculated 

above  instead  of  the  free-stream  value.  The  result  is  a velocity 
component  profile  similar  in  shape  to  the  profile  obtained  on  the  basis 
of  the  free-stream  Mach  number  but  shifted  in  the  upstream  direction. 

In  other  words,  the  point  on  the  xw  traverse  where  the  first  wing 
thickness  effects  are  felt  now  lies  ahead  of  the  point  shown  in  the 
previous  sketch. 

Different  local  Mach  numbers,  M^,  can  be  calculated  for  additional 
traverses  parallel  to  the  x^  traverse  and  located  in  the  same  y plan 
usir\g  the  method  described  in  the  previous  section.  After  recomputing 
the  flow  field  due  to  wing  thickness  along  the  traverses  with  the  appro- 
priate local  Mach  numbers,  M^,  the  points  where  the  first  influence  is 
felt  describe  a locus  which  is  taken  as  the  shock  shape.  A typical 
result  is  indicated  in  the  sketch  below.  Thus,  points  1,  2,  and  3 are 
the  locations  where  wing  thickness  effects  are  first  calculated  on  the 
basis  of  local  Mach  numbers,  M#  , M.  , and  H.  , respectively.  The 


Wing  root  chord 


•Local  chord 


Traversas 
in  vertical 
plane 


Points  where  first 
wing  influence  is 
felt 


Shock 


Mach  line  based 


*< 


Asymptote 
to  shock 


asymptote  of  the  locus  determined  this  way  runs  parallel  to  the  Mach  line 
based  on  the  free-stream  Mach  number. 

The  flow  field  at  a field  point  under  the  wing  including  the  non- 
linear effects  of  the  shock  and  the  rest  of  the  configuration  is  deter- 
mined as  follows.  For  the  given  field  point  the  local  Mach  number,  M^, 
is  calculated  based  on  wing  thickness  only  for  the  traverse  corresponding 
to  the  lateral  coordinates  yw»^  of  the  point  in  accordance  with  equa- 
tion (28) . The  singularities  used  to  flow  model  the  entire  airframe 
subject  to  the  flight  conditions  at  hand  include  the  constant  u-velocity 
panels  on  the  wing,  pylon  (if  present),  and  interference  shell.  Their 
strengths  have  been  calculated  based  on  the  free-stream  Mach  number  from 
the  set  of  simultaneous  equations  generated  by  the  flow  tangency  condi- 
tions, equations  (9),  (10),  and  (11).  The  strengths  of  the  source  panels 

used  to  model  wing  and  pylon  thickness  are  related  to  the  streamwise 
thickness  slopes,  not  Mach  number.  Assuming  that  the  shock  location  is 
not  altered  by  the  effects  of  the  constant  u-velocity  panels,  their 
contribution  to  the  flow  field  is  then  calculated  on  the  basis  of  the 
local  Mach  number,  M^ . To  this  is  added  the  effects  of  the  wing  source 
panels  also  based  on  Finally,  the  contribution  from  the  line 

sources,  sinks  and  line  doublets  computed  using  the  free-stream  Mach 

number,  M„,  are  added. 
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If  a pylon  is  present  and  if  the  field  point  is  positioned  below 
the  pylon,  the  contribution  to  the  flow  field  from  the  constant  u-velocity 
and  source  panels  on  the  pylon  are  calculated  using  the  free-stream  Mach 
number.  This  is  motivated  by  the  evidence  offered  by  the  schlieren 
pictures  taken  during  the  test  program.  They  did  not  show  traces  of 
shocks  from  the  pylon  leading  edge  extending  in  the  downward  direction 
even  though  the  downwash  data  show  significant  effects.  A shock  can 
emanate  from  the  pylon  leading  edge  and  extend  on  either  side  of  the 
pylon  but  its  effects  are  probably  weak  in  the  region  below  the  tip  of 
the  pylon. 

If  flow  field  properties  are  to  be  calculated  at  a series  of  points 
along  an  axial  traverse  under  the  wing,  all  contributions  except  those 
due  to  the  pylon  are  calculated  using  the  method  described  above.  The 
resulting  calculated  flow  field  may  display  a certain  amount  of  scatter 
which  is  due  to  the  discrete  paneling  method  used.  This  behavior  is 
especially  pronounced  when  the  wing  has  a supersonic  leading  edge.  In 
this  case,  the  major  contribution  to  the  flow  field  under  the  wing  is 
associated  with  the  panels  directly  above  the  traverse.  Increasing  the 
number  of  panels  in  a chordwise  row  reduces  the  amount  of  scatter  some- 
what but  increases  computer  storage  and  time  requirements.  The  axial 
loading  distributions  on  a store  calculated  using  velocity  profiles  with 
scatter  are  characterized  by  rapid  oscillations  with  large  magnitude. 

This  is  caused  by  the  fact  that  the  present  loading  calculation  method 
is  based  on  slender-body  theory  and  requires  rates  of  change  with  axial 
distance  of  the  upwash  and  sidewash  velocity  components.  These  rates  of 
change  can  be  very  large  if  the  calculated  axial  profiles  of  the  velocity 
components  exhibit  scatter  and  can  alternate  sign.  In  order  to  alleviate 
this  problem,  the  calculated  flow  field  is  conditioned  in  the  manner 
described  next. 

Over  a finite  length  within  the  region  of  influence  associated  with 
wing  thickness,  a smoothing  procedure  is  applied  to  the  calculated  flow 
field.  First,  the'  finite  length  is  determined.  At  points  along  the 
traverse  under  consideration,  the  effects  of  the  source  panels  on  the 
wing  are  computed  on  the  basis  of  the  local  Mach  number,  M,,  given  by 
equation  (28) . The  desired  length  or  smoothing  range  is  bracketed  by 
the  point  immediately  behind  the  first  influence  and  the  point  immediately 
in  front  of  the  expansion  cone  associated  with  the  trailing  edge.  Over 
the  range,  the  velocities  due  to  all  aircraft  components  except  the  pylon 
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are  calculated.  This  is  the  velocity  field  to  be  smoothed.  The 
smoothing  function  is  a cubic  polynomial  and  the  coefficients  are 
determined  by  an  application  of  the  least-squares  method. 

Ahead  of  the  smoothing  range  (in  fx'ont  of  the  predicted  3hock)  and 
at  points  behind,  the  flow  field  is  calculated  in  the  usual  manner  as 
the  sum  of  the  effects  induced  by  all  singularities.  Inside  the  range, 
the  function  representing  the  smoothed  velocity  profile  accounts  for 
the  effects  induced  by  the  constant  u-velocity  on  the  wing  and  inter- 
ference shell,  source  panels  on  the  wing,  and  the  line  sources,  sinks 
and  line  doublets  modeling  the  fuselage.  Pylon  effects  are  then  cal- 
culated separately  with  the  free-stream  Mach  number  and  added  to  the 
smoothed  profile. 

The  effect  of  the  smoothing  procedure  is  shown  in  figure  5(a) 
which  shows  the  upwash  at  the  one-third  seraispan  and  1.47  inches  below 
the  wing  of  the  configuration  shown  in  figure  6,  The  upwash  is  calcu- 
lated for  5°  angle  of  attack  and  Mach  number  equal  to  2.0  at  points 
along  the  centerline  of  the  store  0.1  inch  below  the  attached  condition. 
The  nose  of  the  store,  x «■  0,  is  17.65  inches  behind  the  fuselage  nose. 

o 

Behind  the  shock,  the  unsmoothed  upwash  (solid  symbols)  displays  a 

considerable  amount  of  scatter.  The  cubic  polynomial  is  applied  to  the 

smoothing  range  which  in  this  case  starts  at  x - 0.3  and  ends  at 

s s 

some  distance  behind  the  store.  The  smoothing  function  produces  the 

solid  line  representing  the  smoothed  upwash  over  the  range  of  interest. 

Figure  5(b)  shows  the  corresponding  calculated  normal-force  distributions. 

The  solid  symbols  represent  the  results  calculated  using  the  upwash 

indicated  in  figure  5(a)  by  the  solid  symbols.  The  result  obtained  with 

the  smoothed  upwash  profile  is  given  by  the  solid  line.  The  former  is 

seen  to  oscillate  wildly.  The  latter  generates  only  one  pronounced  peak 

at  x /£  - 0,275  and  removes  the  oscillations  towards  the  aft  end  of 

s s 

the  store. 

4.  STORE  FORCE,  MOMENT,  AND  TRAJECTORY  CALCULATION  METHODS 

The  methods  and  equations  used  to  calculate  the  aerodynamic  forces 
and  moments  acting  on  the  separated  store  at  each  point  in  its  trajec- 
tory are  unchanged  from  those  presented  in  section  5 of  reference  2. 

The  only  changes  to  the  discussion  presented  there  are  with  regard  to 
the  calculation  of  the  perturbation  velocity  field  and  the  flow  separation 


30 


I 


location.  Also,  all  references  made  in  that-  discussion  to  incompressible 
and  compressible  velocities  or  planes  should  be  disregarded  in  the 
application  of  the  loading  calculation  methods  to  this  work. 

In  order  to  calculate  forces  and  moments  acting  on  a store  during 
its  trajectory,  the  nonuniform  velocity  field  in  which  the  store  is 
immersed  must  be  determined  at  each  point  in  time.  The  velocity  distri- 
bution along  the  body  axis  is  required  for  the  purpose  of  calculating 
forces  and  moments  acting  on  the  body.  The  method  used  to  calculate  the 
flow  field  is  described  in  the  preceding  section.  In  calculating  the 
velocity  field  for  the  case  of  a store  released  from  under  the  fuselage, 
the  procedure  used  is  the  same  as  for  the  under-wing  case.  All  airframe 
components  are  accounted  for  in  contrast  with  the  corresponding  procedure 
described  in  section  5 of  reference  2.  If  an  empennage  is  present,  the 
velocities  also  need  to  be  calculated  at  points  on  the  tail  fin  surfaces 
for  the  empennage  force  and  moment  calculation. 

In  this  supersonic  store  trajectory  program,  the  store  force  and 
moment  methods  are  still  based  on  slcmder-body  theory.  As  will  be  seen 
later,  slender-body  theory  is  not  sufficiently  accurate  for  calculating 
loading  distributions  on  stores  in  cupersonic  nonunifom  flow  fields 
under  the  wing  even  though  the  overall  forces  and  moments  are  predicted 
fairly  well.  The  loading  calculation  methods  presently  implemented  in 
the  computer  program  can  be  upgraded  by  means  of  a three-dimensional 
calculation  method  Such  a refined  method  was  developed  during  the  sub- 
sonic work  and  improved  certain  subsonic  results  appreciably.  It  is 
described  in  detail  in  Appendix  III  of  reference  2. 

The  computer  program  contains  an  input  for  the  axial  location  along 
the  store  where  boundary-layer  separation  occurs.  For  the  subsonic  work, 
it  is  based  on  empirically  obtained  correlations  performed  on  bodies  of 
revolution  at  low  Mach  numbers.  As  such,  equation  (50)  in  reference  2 
is  no  longer  valid  for  the  supersonic  speed  regime  of  this  work  and 
should  be  updated  as  the  program  is  developed  further. 

The  method  used  to  compute  empennage  forces  is  unchanged.  It  is 
based  on  slender-body  theory  with  a correction  for  aspect  ratio  and 
includes  an  application  of  reverse-flow  theorems.  If  a shock  inter- 
sects the  empennage,  strong  chordwise  variations  in  the  flow  field  can 
occur.  The  present  method  may  need  to  be  upgraded  to  a three-dimensional 
lifting-surface  method  to  handle  such  a situation.  In  most  cases,  the 
empennage  is  not  subjected  to  such  a flow  field-. 
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The  trajectory  is  calculated  using  the  six-degree-of-f reedom 
equations  of  motion  used  in  the  subsonic  work.  The  underlying  theory 
is  given  in  full  detail  in  section  6 and  Appendix  II  of  reference  2. 

5.  COMPARISONS  WITH  EXPERIMENTAL  DATA 

This  section  will  present  comparisons  between  results  predicted  with 
the  present  computer  program  and  experimental  data.  The  comparisons 
serve  to  assess  the  accuracy  of  the  supersonic  flow  models  as  well  as 
the  flow  field  and  loading  calculation  methods  described  in  sections  3 
and  4 of  this  report. 

Comparisons  are  made  with  flow  field  data,  store  body  load  distri- 
butions, and  store  body  forces  and  moments.  The  aircraft  configuration 
consists  of  a wing-fuselage  combination  with  and  without  a pylon  attached. 

5.1  Wind-Tunnel  Model  Description 

For  the  purpose  of  providing  systematic  data  with  which  to  assess 
the  accuracy  of  the  present  theory,  a wind-tunnel  test  program  was  con- 
ducted in  the  supersonic  Tunnel  A at  the  Arnold  Engineering  Development 
Center.  The  basic  model  used  in  the  supersonic  tests  was  also  used  in 
the  subsonic  work  reported  in  references  1,  2,  and  5.  This  model  is 
shown  iri  figure  6(a)  . For  some  of  the  teat  results  reported  here  a 
£>ylon  was  attached  to  either  the  wing  or  fuselage.  Figure  6(b)  shows 
the  double-wedge  shape  of  the  pylons.  One  pylon  could  be  attached  to 
the  lower  wing  surface  at  the  one-third  or  two-third  semispan  position 
and  another  could  be  attached  to  the  bottom  of  the  fuselage.  On  the  wing, 
the  pylon  centerline  is  at  40.0  percent  of  the  wing  chord.  Below  the 
fuselage,  the  pylon  leading  edge  is  18.1  inches  behind  the  fuselage  nose. 

The  store  used  in  the  test  program  is  shown  in  figure  6(c) . The 
version  of  the  store  used  in  the  majority  of  the  tests  was  without  the 
cruciform  tail  fins  shown  in  the  figure.  The  pressure  distribution 
model  does  not  have  cruciform  tail  fins  but  is  instrumented  with  19 
pressure  orifices  equally  spaced  in  a meridian  plane  over  the  length 
of  the  store  body.  By  rolling  the  model  through  360°,  a complete  pres- 
sure distribution  was  obtained.  The  load  distributions  were  then 
determined  by  a numerical  integration  of  the  pressures. 

Some  flow  field  comparisons  will  be  shown  which  are  associated  with 
a sharp-edged  trapezoidal  wing  attached  to  an  axisymmetric  fuselage 
above  the  raidwing  position.  This  model  is  described  in  reference  13. 
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5.2  Shock  Shapes 


In  the  section  concerned  with  the  flow  field  calculation  method, 
it  is  mentioned  that  linear  theory  by  itself  does  not  account  for  shock 
strength,  location,  and  curvature.  However,  in  the  same  section  a method 
is  described  which  can  be  viewed  as  partially  accounting  for  these  shock 
phenomena  in  calculating  supersonic  flow  fields  under  the  wing.  As  a 
first  test  of  the  method,  the  predicted  locus  of  points  where  affects 
of  the  wing  are  first  felt  (on  the  basis  of  the  calculated  local  Mach 
number,  M^)  can  be  compared  with  the  shape  deduced  from  the  experimental 
data. 


Comparisons  between  predicted  shock  shapes  and  the  shock  shapes 
obtained  from  the  flow  field  measurements  are  shown  in  figure  7 for  free- 
stream  Mach  numbers  of  1.5,  2.0,  and  2.5.  Data  and  predicted  results  are 
plotted  for  the  one-third  semispan  location  under  the  wing.  A Mach  .line 
is  shown  which  is  based  on  the  free-stream  Mach  nuia1  jr  and  drawn  from  the 
leading  edge  of  the  local  wing  chord.  The  experimental  data  and  the  pre- 
dicted shock  shapes  are  seen  to  be  displaced  forward  from  the  Mach  line. 
All  predictions  generated  by  the  present  computer  program  are  designated 
present  method.  For  all  Mach  numbers  the  left  wing  is  covered  with  a 
constant  u-velocity  panel  layout  consisting  of  8 chordwise  rows  with 
12  panels  in  each  row.  In  addition,  there  are  8 chordwise  rows  with 
20  source  panels  in  each  row  to  account  for  thickness.  The  maximum  values 
used  for  the  thickness  slope  are  determined  as  described  in  section  3.4.1, 


and  the  results  are  tan  0 - 0.209,  0.084,  and  0.186  for  M - 1-5. 

2.0,  and  2.5,  respectively.  On  the  fuselage,  the  interference  shell  is 
laid  out  over  the  length  of  t ha  wing-fuselage  junction  and  is  covered 
over  its  left  half  with  12  half  rings  each  with  4 constant  u-velocity 
panels  above  and  below  the  wing-fuselage  junction.  The  fuselage  is 
modeled  by  53  line  sources,  sinks  and  line  doublets.  The  coordinates 


(x„. z„)  shown  in  figure  7 correspond  to  the  fuselage  coordinate  system 

D H 

of  figure  1 with  the  origin  at  the  nose  tip. 


The  experimentally  obtained  {joints  correspond  to  the  axial  location 
of  the  data  point  immediately  after  the  first  peak  in  the  measured  flow 
fields.  Some  of  the  flow  fields  are  shown  in  the  next  section.  Two 
seta  of  points  are  plotted,  one  set  for  zero  angle  of  attack  and  the 
other  for  5°.  It  is  seen  that  the  shock  location  is  not  influenced 
greatly  by  angle  of  attack  effects.  This  fact  is  used  as  an  assumption 
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in  the  prediction  method  for  calculating  flow  fields  under  the  wing. 

Note  that  the  wing  features  a blunt  leading  edge  for  all  three  cases 
considered  here. 

Figure  7(a)  for  ■ 1.5  corresponds  to  an  almost  sonic  wing 
leading  edge  (component  of  mm  normal  to  the  edge  is  about  1) . The 
shock  appears  to  stand  ahead  of  the  local  wing  chord  shown  at  the  top 
of  the  figure.  In  fact,  the  shock  originates  at  the  leading  edge  of 
the  wing-fuselage  junction  for  this  case.  A slight  dependence  on  angle 
of  attack  is  indicated.  The  predicted  shock  shape  indicated  by  the 
solid  line  matches  the  shock  shape  obtained  from  experiment  well. 

The  case  for  Mw  •<  2 is  shown  in  figure  7(b)  and  corresponds  to  a 
supersonic  wing  leading  edge  (component  of  normal  to  the  edge  is 

greater  than  1) . Except  for  the  data  point  at  the  lowest  position 
(zfi  = 4.37  inch« s) , the  experimental  shock  shapes  for  af  - 0°  and  5° 
coincide.  The  present  method  predicts  a shock  shape  which  lies  slightly 
aft  of  the  experimental  one.  The  agreement  is  considered  good.  An 
additional  theoretical  curve  is  given  by  the  dashed  line.  It  represents 
a hybrid  shock-shape  calculation  approach  described  in  Appendix  111.  In 
short,  this  calculation  combines  a simple  concept  and  curve  fitting 
procedure  of  Love  described  in  reference  14  with  input  for  shock  standoff 
distance  generated  by  the  numerical  approach  of  Rizzi  discussed  in  ref- 
erence 15.  The  hybrid  theory  indicates  a shock  shape  slightly  ahead  of 
the  experimental  one  and  definitely  ahead  of  the  shape  predicted  by  the 
present  method.  As  an  approximate  indication  of  the  difference  between 
pure  linear  theory  and  the  present  method,  the  Mach  line  based  on  M^ 
and  with  its  origin  on  the  wing  leading  edge  is  shown. 

Figure  7(c)  3hows  the  same  set  of  results  for  *»■  2.5.  The  leading 
edge  of  the  wing  is  supersonic  again,  and  the  shock  has  greater  sweepback 
than  that  for  the  Mw  ■ 2 results.  This  trend  is  also  indicated  by  the 
prediction  from  the  present  method.  For  this  flow  condition  the  present 
method  is  represented  by  the  solid  line  on  the  figure  faired  through 
five  calculation  points,  the  highest  of  which  is  at  zB  » 1.37  inches. 

As  it  stands,  agreement  with  experiment  is  good  and  the  result  from  the 
hybrid  theory  now  lie3  behind  the  data  and  the  present  method.  The 
shock  shape  indicated  by  experiment  shows  a small  effect  due  to  angle  of 
attack. 


5.3  Flow  Fields 


The  capability  of  the  present  method  for  predicting  shock  shapes  in 
the  vertical  plane  is  discussed  above.  It  remains  to  be  seen  how  well 
the  present  method  can  predict  the  flow  field  accounting  for  the  non- 
linear shock  behavior.  To  assess  this,  extensive  comparisons  with 
experimental  flow  field  data  from  the  test  program  will  be  shown.  The 
wind-tunnel  model  is  the  wing-fuselage  combination  of  figure  6(a).  The 
effects  of  the  pylon,  shown  in  figure  6(b),  will  also  be  discussed. 
Comparison  with  data  from  a configuration  equipped  with  a sharp  leading- 
edge  wing  is  also  shown.  All  predictions  shown  are  based  on  the  same 
panel  layout  used  for  the  determination  of  the  shock  shape  as  discussed 
previously  including  the  specification  for  maximum  thickness  slope. 

When  a pylon  is  attached,  a layout  of  2 chordwise  rows  with  4 constant 
u-velocity  panels  in  each  row  is  used.  To  model  thickness,  the  source 
panel  layout  on  the  pylon  consists  of  2 chordwise  rows  with  9 panels 
in  each  row.  The  thickness  slopes  are  either  constant  or  zero,  if  the 
chordwise  length  of  the  source  panels  could  be  varied,  only  3 source 
panels  would  be  required  to  model  pylon  thickness.  At  present,  the 
computer  program  allows  only  for  equally  spaced  panels  in  the  chordwise 
direction.  Along  the  span,  the  number  of  source  panels  must  equal  the 
number  of  constant  u-velocity  panels.  For  thi3  type  of  thickness  distri- 
bution it  would  be  advantageous  to  be  able  to  specify  the  source  panels 
layout  independently. 

5.3.1  Results  for  Mach  number  1.5 

The  effects  of  vertical  distance  from  the  wing  are  shown  in  fig- 
ures 8 through  10  for  zero  angle  of  attack  and  free-stream  Mach  number 
equal  to  1.5.  The  coordinates  (xB,yB,zQ)  shown  on  the  sets  of  figures 
pertain  to  the  fuselage  coordinate  system  with  the  origin  at  the  nose 
tip  as  shown  in  figure  1.  The  flow  field  is  characterized  by  the  fol- 
lowing set  of  velocity  components.  Backwash  u i3  positive  i,u  the 
negative  xB  direction.  Sidewash  v is  positive  in  the  inboard 
direction  under  the  left  wing  and  upwash  w is  positive  upwards.  The 
vertical  location  of  the  traverse  and  the  local  wing  chord  are  indicated 
at  the  top  in  each  figure.  Wing  thickness  is  not  drawn  to  scale.  For 
all  vertical  locations,  the  leading-edge  shock  and  trailing-edge  Mach 
cone  associated  with  the  wing  are  indicated  by  the  sudden  jump  in 
magnitude  of  the  measured  backwash,  sidewash,  and  upwash.  For 
zB  ~ 4.37  inches  (fig.  10) , the  trailing-edge  Mach  cone  effect  is  beyond 
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the  range  of  the  traverse  shown  here.  The  present  method  predicts  the 
jumps  well  for  all  vertical  locations  although  the  wing  trailing-edge 
effect  is  shifted  forward  slightly.  It  is  interesting  to  note  that  the 
upwash  and  to  a lesser  extent  the  backwash  are  still  strongly  felt  at 
the  lowest  (zB  ■*  4.37  inches)  position,  figures  10(a)  and  10(c),  whereas 
the  sidewash  is  almost  completely  attenuated  in  figure  10(b).  This 
confirms  the  fact  that  the  main  effect  of  the  3hock  wave  is  to  impart  a 
strong  downward  change  in  momentum  of  the  flow.  The  present  method  pre- 
dicts the  characteristics  of  the  flow  field  very  well  except  immediately 
behind  the  shock.  The  theory  tends  to  ovarpredict  the  magnitudes  of  the 
downwash  immediately  behind  the  forward  jump  or  shock  wave.  According 
to  the  observations  reported  in  reference  16,  shock  waves  tend  to  soften 
sharp  peaks  in  static  pressure  signatures  when  theoretical  discontinuities 
in  the  slopes  are  indicated. 

Only  the  sidewash  and  upwash  velocity  profiles  are  used  in  the 
loading  calculation  method  based  on  slender-body  theory  mentioned  in 
section  4.  Therefore,  in  the  following  comparisons  the  backwash  will 
not  be  shown  and  discussed.  Generally,  the  agreement  between  theory 
and  experiment  in  backwash  is  similar  to  the  agreement  indicated  for 
upwash. 

Figures  11,  12,  and  13  show  the  variation  of  the  sidewash  and  upwash 
with  vertical  distance  for  5°  angle  of  attack  and  a free-stream  Mach 
number  1.5.  The  figures  show  the  3ame  trends  as  shown  in  figures  8,  9, 
and  10.  As  observed  for  the  zero  angle  of  attack  case,  the  measured 
sidewash  attenuates  more  rapidly  than  the  measured  upwash.  For 
Zjj  w 4.37  inches,  the  latter  has  a jump  in  magnitude  equal  to  about 
70  percent  of  the  jump  for  za  - 1,37  inches.  This  behavior  is  well 
predicted.  For  the  two  higher  locations  (zQ  «■  1.37  and  2.12  inches), 
the  present  method  calculates  the  upwash  better  than  the  sidewash.  The 
location  of  the  shock  wave  is  predicted  about  1/2  inch  behind  the  location 
indicated  by  the  data  for  zn  - 2.12  inches.  This  fact  is  also  indicated 

O 

on  figure  7 (a) . 

It  is  now  of  interest  to  study  the  capability  of  the  present  method 
to  account  for  the  addition  of  the  pylon  shown  on  figure  6(b)  to  the 
basic  wing-fuselage  configuration  of  figure  6(a).  The  results  are  shown 
in  figure  14  for  0°  angle  of  attack  and  in  figure  15  for  5°.  At  the  top 
of  these  figures,  the  wing  local  chord  and  pylon  are  indicated.  Also 
shown  are  two  Mach  lines,  based  on  the  free-stream  Mach  number,  extending 
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from  the  pylon  leading  and  trailing  edge.  These  lines  serve  to  indicate 
the  region  of  pylon  influence,  especially  the  onset  of  pylon  effects. 

If  the  assumption  under  which  pylon  effects  arc  accounted  for  is  correct, 
the  pylon  effects  should  be  first  felt  at  the  intersection  of  the  forward 
Mach  line  and  the  traverse.  The  assumption  refers  to  the  region  of  influence 
of  the  pylon  being  determined  by  the  frae-stream  Mach  number.  In  both 
cases,  the  upwash  is  affected  strongly  by  the  addition  of  the  sharp-edged 
pylon.  The  present  method  predicts  the  levels  well  although  the  pre- 
dicted onset  of  pylon  effect  is  slightly  behind  the  experimental  results 
for  the  5°  angle  of  attack  case.  It  is  interesting  to  compare  these 
data  with  that  obtained  during  the  subsonic  work  with  a pylon  with 
shorter  chord  and  blunt  leading  and  trailing  edges.  Such  data  are  given 
in  figure  14  of  reference  2.  In  comparison,  the  subsonic  pylon  effects 
are  similar  in  nature  but  much  more  subdued  notwithstanding  the  blunt 
edges.  In  general,  for  supersonic  flow  the  flow  features  tend  to  be 
more  peaked. 

The  pylon  trailing  edge  is  at  - -22.2  inches  and  for  axial 
locations  beyond  - -24.5  inches,  the  predicted  effect  of  the  pylon 
on  sidewash  does  not  occur,  and  the  effect  of  the  pylon  washes  out 
experimentally.  This  result  is  believed  due  to  flow  separation  on  the 
aft  portion  of  the  pylon.  It  was  found  that  in  this  region  the  dominant 
contribution  to  the  predicted  results  is  generated  by  the  constant 
u- velocity  panels  on  the  wing  in  the  chordwise  rows  adjacent  to  the 
pylon.  In  particular,  the  aft  panels  were  heavily  influenced  by  the 
pylon  source  panels  modeling  thickness.  Thus,  if  complete  or  partial 
flow  separation  occurred  on  the  back  portion  of  the  pylon,  the  pylon 
thickness  interference  on  the  wing  panels  would  be  reduced  or  possibly 
completely  cancelled.  As  a result,  only  the  wing  source  panels  would 
contribute  and  the  prediction  would  fall  on  the  pylon-off  (dashed)  line. 
The  downwash  does  not  appear  to  be  influenced  this  way.  The  overall 
agreement  between  predicted  and  measured  pylon  effect  on  upwash  i3  good. 

For  af  ■*  0°,  the  measured  rear  peak  in  upwash  shown  in  figure  14(b) 
lies  slightly  aft  of  the  predicted  peak  while  for  “ 5°  the  reverse 

is  true.  In  this  region  of  the  flow  field,  the  flow  is  accelerated  for 
the  zero  angle  of  attack  case  causing  the  local  Mach  number  to  be  higher 
than  calculated  behind  the  shock  and  possibly  even  higher  than  the 

free-stream  value,  Mm.  Thus,  the  direct  pylon  effects  may  travel  along 
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Mach  lints  swept  hack  more  than  ths  ones  shown  at  the  top  of  the  figure. 

The  present  method  computes  the  pylon  influence  on  the  basis  of  free- 
stream  Mach  number,  M^.  For  af  ■ 5°,  the  local  Mach  number  in  the 
region  where  the  rear  peak  occurs  may  not  be  as  high  as  Mw  since  the 
flow  is  compressed  below  the  lifting  wing.  Therefore,  the  pylon  effects 
may  travel  along  Mach  lines  with  less  sweep  than  the  lines  shown  at  the 
top  of  figure  15(b) . To  account  for  this  nonlinearity  requires  a further 
correction  to  the  one  presently  employed  for  the  wing  leading-edge  shock 
but  may  not  be  worth  the  additional  effort  it  entails.  Increasing  the 
number  of  constant  u-velocity  panels  on  the  pylon  may  influence  the 
results.  All  results  presented  here  are  based  on  a total  of  200  constant 
u-velocity  panels  and  the  present  computer  program  is  limited  to  200 
panels. 

The  effects  of  the  pylon  on  the  downwash  under  the  fuselage  are 
shown  in  figure  16  for  0°  and  5°  angle  of  attack.  Compared  to  the 
experimental  data  under  the  wing,  these  data  exhibit  an  appreciable 
amount  of  scatter  especially  with  the  pylon  removed.  In  part,  this  is 
suspected  to  be  due  to  a flat  surface  to  which  the  pylon  can  be  mounted 
which  is  recessed  slightly  into  the  curved  underside  of  the  fuselage. 

For  a«j  <■  0,  agreement  between  present  method  and  experiment  ia  consid- 
ered good.  The  measured  pylon  effects  are  of  the  same  magnitude  as  for 
the  under-wing  case  and  are  predicted  well  by  the  present  method.  In 
figure  16(b),  the  present  method  predicts  the  minimum  value  of  the 
pylon-on  upwash  profile  in  part  only  but  in  general  the  agreement  ia 
considered  good. 

5.3.2  Results  for  traperoidal  wing-fuaelaga  combination  for  Mach 
number  1.61 

For  the  purpose  of  testing  the  capability  of  the  present  method  for 
predicting  flow  fields  under  a sharp-edged  wing,  experimental  data  ware 
taken  from  reference  13.  The  configuration  with  pertinent  specifications 

i i 

w 1 is  shown  in  the  following  sketch  taken  from  that  reference.  All  dimensions 
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0.69  2.75 


Trapezoidal  Wing-Fuselage  Configuration 
Aspect  ratio  3.0 
Taper  ratio  0.25 

Section  4-percent  circular  arc 

Area  1.00  square  foot 


shown  are  in  inches.  The  plane  of  the  wing  is  about  one-half  of  the 
fuselage  radius  above  the  fuselage  midplane. 

For  the  free-stream  Mach  number  of  1.61,  the  wing  leading  edge  is 
supersonic.  The  maximum  streanrwise  thickness  slope  occurs  at  the 
leading  and  trailing  edges  and  its  magnitude  is  0.076.  This  slope 
corresponds  to  an  angle  of  4.36°  measured  in  the  streamwisa  direction. 

In  the  direction  normal  to  the  wing  leading  edge,  equation  (21)  gives 
an  angle  of  about  5.1°.  The  same  equation  ulso  gives  the  Mach  number 
normal  to  the  leading  edge,  M^,  to  be  equal  to  1.38.  After  consulting 
figure  4 in  reference  11  or  the  more  precise  chart  2 in  the  same  ref- 
erence, the  maximum  wedge  semi-vertex  angle  for  an  attached  shock  is 
7.8°.  Thus,  the  shock  is  attached  to  the  wing  leading  edge. 

Figures  17  and  18  are  comparisons  between  the  results  of  the 
present  method  and  experimental  data  drawn  from  reference  13.  The 
results  are  expressed  as  local  upwash  angle,  and  sidewash  angle,  cr^. 

In  terras  of  the  perturbation  velocities  uw,vw,ww  in  the  wing  coordinate 
system  of  figure  1,  the  angles  are  expressed  as  follows. 
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At  the  top  of  both  figures,  the  local  wing  chord  and  vertical  location 
of  the  traverse  are  indicated.  Wing  thickness  is  not  drawn  to  scale. 
Mach  lines  based  on  M^  are  also  shown  to  indicate  approximately  the 
linear  region  of  influence. 


The  results  for  zero  angle  of  attack  are  given  in  figure  17. 

Overall,  the  present  method  predicts  the  sidewash  and  upwash  well,  and 
the  systematic  differences  in  the  neighborhood  of  the  leading-edge  shock 
formerly  seen  are  exhibited  here  also.  However,  the  measured  shock 
location  lies  about  1/2  inch  forward  of  the  predicted  shock  location. 

The  reason  for  this  discrepancy  is  presently  not  known.  The  wing  leading 
edge  may  not  be  aerodynamically  sharp;  that  is,  it  may  be  blunt  enough 
to  cause  the  shock  to  be  slightly  detached.  On  the  other  hand,  the 
calculation  of  the  local  Mach  number  in  accordance  with  equation  (28) 
can  be  sensitive  to  the  spanwise  number  of  source  panels  laid  out  on  the 
wing.  The  shock  location  is  determined  by  the  local  Mach  number.  The 
present  method  results  are  based  on  the  same  panel  layout  as  specified 
in  an  earlier  section  concerned  with  shock  shapes.  Consequently,  the 
results  generated  by  the  present  method  are  based  on  8 chordwise  rows 
with  20  source  panels  in  each  row.  A better  local  Mach  number  may  result 
with  more  chordwise  rows,  but  this  has  not  been  ascertained  as  yet. 


Figure  18  shows  the  results  for  ct^  — 4°.  The  trend  of  sidewash 
and  upwash  with  jangle  of  attack  are  predicted  well.  For  example,  the 
measured  sidewash  is 'definitely  increased  in  the  outboard  direction  as 
the  angle  of  attack  is  increased  to  4°.  The  present  method  also  indi- 
cates this  behavior.  Otherwise,  the  same  remarks  made  above  for  the 
zero  angle  .of  attack  case  apply. 
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5.3.3  Results  for  Mach  number  2 


The  effects  of  the  pylon  on  the  flow  field  under  the  wing  are  shown 
in  figures  19  and  20  for  0°  and  5°  angles  of  attack  respectively.  At 
the  top  of  the  figures,  the  local  wing  chord  and  the  vertical  location 
of  the  traverse  are  indicated.  Wing  thickness  is  not  drawn  to  scale. 

The  pylon  is  also  outlined  and  free-streara  Mach  lines  are  drawn  from 
the  leading  and  trailing  edges  of  the  pylon  tip. 

For  both  angles  of  attack,  the  predicted  shock  location  lies  behind 
the  experimentally  deduced  location  by  about  1/2  inch.  With  the  pylon 
removed,  the  upwash  is  predicted  better  than  the  sidewash.  Generally, 
speaking,  the  pylon-off  velocity  profiles  are  predicted  reasonably  well, 
except  for  the  measured  sidewash  and  downwash  values  at  the  leading-edge 
shock . 

In  comparison  with  the  experimental  data  for  Mach  number  1.5  shown 
in  figure  14,  the  effects  of  the  pylon,  shown  in  figure  19,  are  somewhat 
less  pronounced  for  Mach  number  2.  The  backwards  shift  of  the  location 
of  the  shock  is  also  evident.  The  present  method  for  af  ■ 0°  predicts 
these  trends  fairly  well.  An  increase  in  sidewash  due  to  the  pylon  of 
about  1°  is  indicated  by  the  present  method  at  points  behind  x^  ■»  -23 
inches  whereas  the  measurements  only  show  a fractional  increase. 

Comparing  the  measured  upwash  profile  of  figure  14(b)  for  Mach  number 
1.5  to  the  profile  for  Mach  number  2 in  figure  19(b)  indicates  very 
little  shift  in  the  location  where  the  pylon  influence  is  first  felt. 

The  theory  indicates  a similar  result.  Otherwise,  the  remarks  made  in 
connection  with  the  Mach  number  1.5  comparison  applies  here  as  well. 

When  the  angle  of  attack  is  changed  to  5°,  the  pylon  exhibits  loading 
due  to  outboard  sidewash  induced  by  the  lifting  wing.  Below  the  pylon, 
the  sidewash  must  then  accelerate  around  the  pylon  tip  in  the  outboard 
direction.  This  increment  in  sidewash  with  angle  of  attack  is  indicated 
by  both  experimental  data  and  present  method  in  figure  20(a).  This 
increment  is  larger  than  the  increment  shown  for  Mach  number  1.5  in 
figure  15(a).  The  trends  shown  by  the  upwash  are  predicted  fairly  well. 
The  data  with  the  pylon  on  do  not  go  back  far  enough  to  show  if  the  pylon 
effect  washes  out  behind  its  trailing-edge  location  as  a result  of 
separation.  Such  separation  effects  would  be  greatly  reduced  at  full- 
scale  Reynolds  number. 
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5.3.4  Results  for  Mach  number  2.5 

This  is  the  highest  Mach  number  for  which  data  were  taken  on  the 
basic  wing-fuselage  model  shown  in  figure  6(a) . Only  the  pylon-off  data 
will  be  discussed.  The  sidewash  and  upwash  comparisons  are  shown  in 
figure  21  for  zero  angle  of  attack,  QSf  “ 0°,  and  in  figure  22  for 
af  - 5°.  At  the  top  of  the  figure,  the  local  wing  chord  and  the  vertical 
location  of  the  traverse  are  shown.  The  indicated  wing  thickness  is  not 
drawn  to  scale.  Mach  lines  based  on  free-streara  Mach  number,  M , are 
drawn  from  the  wing  leading  and  trailing  edges. 

Considering  the  forward  peak  in  sidewash  in  figures  8(b),  19(a), 
and  21(a),  it  is  interesting  to  note  that  the  actual  measured  peak  value 
is  almost  constant  for  Mw  ■ 1.5,  2.0,  and  2.5.  At  the  highest  Mach  number, 
the  measured  and  predicted  axial  location  of  the  peak  are  still  ahead 
of  the  axial  location  of  the  intersection  of  the  Mach  line  based  on  M^ 
and  the  traverse  as  would  be  expected  and  the  shock  wave  is  bent  back 
with  increasing  Mach  number.  The  predictions  generated  by  the  present 
method  are  based  on  the  same  panel  layout  specified  in  the  section 
describing  the  shock  shapes.  The  maximum  thickness  slope  for  this  Mach 
number  is  tan  “ 0.186.  Immediately  following  the  shock,  the 

pre  ant  method  overestimates  the  magnitude  of  the  sidewash  in  a way  which 
tends  to  become  more  pronounced  as  Mw  or  af  increase.  In  terms  of 
upwash,  the  effect  of  higher  Mach  number  is  to  deepen  the  forward  peak 
and  to  shift  back  the  axial  location  of  the  peak  as  can  be  seen  on  fig- 
ures 8(c),  19(b),  and  21(b).  This  is  also  indicated  by  the  present 
method.  For  the  highest  Mach  number,  the  present  method  overpredicta 
the  magnitudes  of  the  upwash  in  the  region  immediately  following  the 
shock,  in  general,  the  velocity  profiles  are  predicted  fairly  well. 

Figure  22  shows  comparisons  for  - 5°.  The  magnitudes  of  side- 
wash  immediately  behind  the  shock  is  overestimated  by  the  present  method. 
This  effect  can  also  be  discerned  for  the  predictions  of  sidewash  in 
figures  15(a)  and  20(a)  for  Mw  = 1.5  and  2,  respectively,  for  the  pylon- 
off  case.  In  fact,  the  measured  sidewash  profiles  show  a systematic 
change  with  M^  for  reasons  which  are  not  known.  The  measured  upwash 
profile  changes  its  shape  behind  the  shock  from  a dish  to  a ramp.  This 
result  is  predicted  by  the  present  method  as  can  be  seen  on  figures  15(b) , 
20(b),  and  22(b).  For  Mm  - 2.5,  agreement  between  the  present  method 
and  measurement  is  not  so  good  as  for  the  lower  Mach  numbers. 
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5.4  Loadings  on  the  Store 

A few  comparisons  between  measurements  and  results  from  the  present 
method  in  terms  of  store  force  distributions  and  overall  forces  and 
moments  will  now  be  shown.  These  comparisons  point  out  the  limitations 
associated  with  determining  the  load  distribution  at  supersonic  speedo 
with  the  present  method  taken  from  the  subsonic  work  of  reference  2 
and  also  suggest  the  means  for  overcoming  these  limitations.  (Neverthe- 
less, the  total  forces  will  be  seen  to  be  good  enough  for  preliminary 
design  purposes.)  A few  pressure  distributions  along  the  length  of  the 
store  body  calculated  by  the  present  theory  are  compared  with  measurements 
for  the  finless  store  in  uniform  flow. 

5.4.1  Force  distributions  along  the  store 

The  results  generated  by  the  present  method  are  based  on  the  panel 
layout  and  other  specifications  listed  in  an  earlier  section  dealing 
with  shock  shapes.  As  formulated  for  use  at  subsonic  speeds,  the  present 
method  calculates  store  loads  on  the  basis  of  slender-body  theory  with  a 
two-dimensional  approach  for  the  buoyancy  forces. 

The  side-force  and  normal -force  distribution  along  the  axis  of  the 
attached  store  are  shown  in  figure  23  for  zero  angle  of  attack,  af  ■ 0°, 
and  for  free-stream  Mach  number  1.5.  These  data  are  obtained  from  the 
pressure  distribution  model  by  an  integration  of  the  pressures.  The 
effect  of  the  pylon  on  the  measured  normal -force  distribution  is  to 
enlarge  and  slightly  shift  back  the  peaks  for  the  pylon-off  case.  The 
aide-force  distribution  is  only  slightly  affected  by  addition  of  the 
pylon.  For  comparison  with  earlier  figures,  the  body  axis  coordinate 
Xg  is  also  shown  in  figure  23. 

in  connection  with  the  side-force  distribution,  the  calculated 
velocity  profile  from  which  the  necessary  gradients  are  computed  is 
almost  identical  to  the  calculated  sidewash  distribution  shown  in  fig- 
ure 14(a) . Apparently,  the  inclusion  of  store  effects  in  the  wing-pylon 
boundary  condition  has  little  effect  on  the  sidewash.  The  store  nose  is 
at  the  location  of  the  3hock  wave  from  the  wing  leading  edge.  The 
comparison  between  the  prediction  and  data  on  the  store  nose  illustrates 
the  shortcoming  of  slender-body  theory  which  is  based  on  small  axial 
gradients  and  as  such  it  cannot  cope  with  the  large  gradients  behind  the 
shock  even  though  they  are  smoothed.  The  sudden  change  in  the  experimental 
side-force  data  at  x /i  approximately  equal  to  0.6  is  definitely  a 
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flow  field  effect  not  accounted  for  by  tbe  present  method  and  is  thought 
to  be  an  extraneous  disturbance  due  to  some  unknown  cause.  It  can  also 
be  seen  in  figure  14(a)  at  Xg  - -22.25  inches.  This  effect  is  not 
caused  by  the  pylon. 

The  normal-force  distribution  shown  in  figure  23(b)  is  strongly 

influenced  by  the  pylon.  From  x /i  - 0.2  to  0.32,  both  pylon-off  and 

pylon-on  experimental  data  show  the  effect  of  the  store  nose  shock 

reflected  from  the  wing.  Beyond  this  range,  the  pylon-on  data  are 

strongly  influenced  by  the  pylon  effects.  For  the  prediction  of  the 

pylon-off  case,  it  was  found  that  the  calculated  but  unsmoothed  velocity 

profile  exhibited  a negative  peak  at  x^  - -20.3  inches  which  was  due  to 

reflection  of  the  store  nose  wave.  Thi3  location  corresponds  to  about 

x /f  ->  0.35  which  is  slightly  behind  the  experimental  negative  peak.  The 
s s 

pylon-off  prediction  for  normal-force  distribution  does  not  show  any 

effect  of  the  reflected  shock  at  x /£  « 0.35.  The  reason  for  this 

s s 

defect  is  that  the  smoothing  procedure,  presently  used  to  condition  the 
velocity  profile  for  use  with  slander-body  theory,  completely  washes  out 
the  effects  of  the  reflected  store  nose  shock.  It  is  interesting  to  note 
that  the  pylon-off  data  resemble  the  pylon-on  data  and  that  the  net  effect 
of  the  peaks  results  in  a small  force  but  appreciable  couple  contribution. 
After  the  negative  peak  at  x /£  ■ 0.32,  the  normal-force  distribution 

rises  to  a positive  peak  and  then  falls  back  on  the  prediction.  The 
store  axis  is  located  1.47  inches  below  the  chordal  plane  of  the  wing. 
Accounting  for  wing  thickness  and  the  store  radius  leaves  a distance  of 
0.78  inch  between  the  under  surface  of  the  wing  and  the  store  body. 
Compared  to  the  local  wing  chord  of  about  7 inches,  this  spacing  makes 
the  store  closely  coupled  with  the  wing  and  multiple  shocks  can  occur 
between  the  wing  and  3tore.  Certain  phenomena  which  can  take  place  in 
this  region  will  now  be  examined. 

Consider  the  following  sketch,  in  accordance  with  reference  17, 
the  parameters  governing  the  reflection  of  shock  waves  from  solid 
boundaries  are  boundary-layer  character  and  thickness,  shock  strength 
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and  the  geometry  of  the  three-dimensional,  shock -hody  configuration. 
Assume  for  now  that  the  shock  from  the  wing  leading  edge  does  not  inter- 
act with  the  store  nose  and  that  the  boundary  layer  on  the  under-wing 
surface  i3  laminar  near  the  leading  edge.  The  boundary  layer  may  then 
separate  at  or  near  the  point  of  impingement  of  the  shock  from  the  store 
nose.  This  shock  is  then  actually  reflected  as  an  expansion  wave.  But 
the  outward  compression  waves  generated  by  the  thickening,  separation, 
and  the  raattachment  of  the  boundary  layer  on  the  wing  tend  to  coalesce 
with  the  expansion  wave  to  form  a compression  shock  marked  1.  Shock  1 
can  now  be  considered  as  the  store  nose  shock  reflected  from  the  wing. 

In  inviscid  flow  terms,  as  this  shock  impinges  on  the  store  body, 
flow  is  induced  into  the  upper  surface  of  the  body.  Since  this  surface 
is  a solid  boundary,  the  required  cancellation  (outward  flow)  brings 
about  another  reflected  shock  2.  As  the  shock  marked  1 intersects  the 
body  on  its  way  to  the  leeward  side,  the  flow  may  leave  the  surface.  The 
cancellation  of  the  outward  component  of  the  flow  on  the  leeward  side 
causes  an  expansion  wave  to  originate  along  the  intersection  line.  It 
can  weaken  shock  1 to  the  point  where  shock  1 is  broken  up.  viscous 
effects  tend  to  distort  the  flow  further  and  in  fact  a very  strong 
boundary-layer  crossflow  occurs  approximately  in  the  direction  of  the 
intersection  line.  The  accumulation  of  boundary  layer  on  the  leeward 
side  regenerates  a compression  shock  some  distance  away  from  the  body 
surface . 

The  crossflow  boundary  layer  has  the  effect  of  relieving  the  over 
pressures  on  the  incident  side  by  allowing  it  to  bleed  off  to  the  leeward 
side.  In  the  absence  of  other  shocks  impinging  on  the  body,  the  over 
pressures  fall  back  to  the  body-alone  pressures  within  a short  distance 
downstream  from  the  inpingement  location.  Thus,  the  reflected  shock 
definitely  causes  a down  force  on  the  body  locally.  Downstream,  the 
local  down  force  should  disappear  and  the  loading  distribution  should 
be  determined  by  the  flow  field  as  seen  by  the  store.  The  pylon-off 
data  shown  in  figure  23(b)  exhibit  a definite  positive  peak,  however. 

It  is  possible  that  the  increasing  upwash,  shown  in  figure  14(b) , 
together  with  the  downward  crossflow  effects  discussed  above,  can 
generate  a complete  reversal  in  the  pressure  field  across  the  body. 
Instead  of  this  phenomenon  or  in  addition  to  it,  shock  2 can  reflect 
off  the  boundary  layer  on  the  wing.  If  the  boundary  layer  is  separated 
due  to  the  action  of  the  store  nose  shock,  shock  2 can  be  reflected  from 
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the  wing  as  an  expansion  wave  towards  the  store  body.  On  the  incident 
side  of  the  store  body,  under  pressures  are  induced  due  to  impingement 
of  the  expansion  wave  and  a positive  normal -force  peak  would  result. 

The  effects  of  the  multiple  shock  system  eventually  die  out  over  the 
aft  portion  of  the  store.  This  trend  is  indicated  by  the  experimental 
data  and  the  present  method  predicts  the  level  of  the  normal-force 
distribution  well.  Thus,  over  the  part  of  the  store  influenced  by  the 
shocks,  the  forces  acting  on  the  store  body  are  dominated  by  the  shocks. 

The  pylon-on  data  appear  to  be  in  phase  with  the  pylon-off  (multiple 
shock)  data.  The  forward  negative  peak  with  the  pylon  on  at  xa/^a  " 0.37 
is  believed  to  be  caused  by  the  pylon  and  is  predicted  well.  The  pylon-on 
prediction  does  not  contain  effects  of  the  reflected  store  nose  Bhock 
because  of  the  present  velocity  smoothing  scheme  as  mentioned  earlier. 

The  pronounced  negative  peak  predicted  by  the  present  method  for  the 
pylon-on  condition  over  the  aft  position  of  the  store  is  due  to  the 
calculated  pylon  induced  flow  field.  The  calculated  upwash  profile  shown 
in  figure  14(b)  shows  a sharp  positive  peak  at  x^  - -22. 5 inches  which 
lies  slightly  ahead  of  the  experimental  one.  Beyond  this  location,  the 
prediction  drops  sharply  thereby  generating  a very  large  negative  gradient. 
Again,  the  magnitude  of  this  axial  gradient  is  not  compatible  with  the 
assumptions  of  alender-body  theory.  Besides,  the  experimental  flow  field 
showe  a less  drastic  drop  off  in  upwash.  Thus,  the  slender-body  load 
prediction  is  also  affected  by  the  inaccuracy  in  the  predicted  upwash 
distribution. 

The  comparisons  exhibited  by  this  figure  suggest  that  better  account 
must  be  taken  of  the  effect  of  the  store  nose  shock  reflection  from  the 
wing  under  surface.  For  stores  close  to  the  wing  where  a number  of 
reflections  can  take  place,  it  is  difficult  to  obtain  the  necessary 
resolution  to  account  for  the  phenomena  without  using  very  large  numbers 
of  wing  constant  u-veiocity  panels  on  the  wing  near  the  store.  However, 
by  using  an  image  store  technique,  it  is  possible  to  account  for  wing 
reflection  effects  very  conveniently  for  those  cases  where  the  store 
nose  shock  is  reflected  from  the  under-wing  surface  and  does  not  pass 
upward  in  front  of  the  leading  edge.  Furthermore,  this  technique  can 
be  combined  with  a three-dimensional  store  loading  calculation  method 
instead  of  the  slender-body  method  and  can  be  carried  out  independently 
of  solving  the  airframe  modeling  problem.  Accordingly,  it  can  probably 
be  carried  out  at  various  points  in  the  store  trajectory  without  using 
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excessive  computer  time.  The  refined,  three-dimensional  store  loading 
method  was  developed  for  subsonic  speeds  in  reference  2.  In  general, 
the  same  type  of  behavior  of  both  prediction  and  measurement  exist  for 
the  other  Mach  numbers. 

5.4.2  Store  forces  and  moments 

The  forces  and  moments  acting  on  the  finlesa  3tore  of  figure  6(c) 
are  shown  in  figure  24  for  5°  angle  of  attack  and  Mach  number  1.5. 

Vertical  location  Az  ■ 0 corresponds  to  the  attached  store  position. 

The  results  shown  are  obtained  with  the  present  method  which  is  based 
on  slender-body  theory.  Note  that  the  measured  side  force  shows  a 
pylon  effect  for  the  first  inch  of  vertical  distance,  but  the  present 
method  does  not  predict  this  trend.  Predicted  side  force  and  yawing 
moment  for  the  pylon-off  case  indicate  the  trends  shown  by  the  experi- 
mental data.  The  normal  force  and  pitching  moment  show  larger  influence 
from  the  pylon  than  do  the  side  force  and  yawing  moment  which  is  also 
indicated  by  the  present  method  for  distances  larger  than  1.5  inches 
below  the  attached  position.  On  an  overall  basis,  experiment  and 
prediction  show  the  same  trends. 

The  behavior  for  Mach  number  2 is  shown  in  figure  25.  Agreement  be- 
tween the  present  method  and  experiment  is  similar  to  the  Mw  ■ 1.5  case 
and  actually  slightly  better  for  the  normal  force  and  pitching  moment 
except  just  below  the  attached  position.  The  prediction  at  Az  - 3 inches 
for  the  side  force  and  yawing  moment  is  excessive.  The  smoothed  sidewash 
velocities  calculated  at  the  rear  of  the  store  resulted  in  large  axial 
changes  which  in  conjunction  with  the  slender-body  load  calculation  method 
generated  unrealistically  high  side  force  and  yawing  moment.  Otherwise, 
the  present  method  indicates  the  measured  trends. 

5.4.3  Finlesa  store  in  uniform  flow 

Tha  purpose  of  the  following  discussion  is  to  compare  loading 
distributions  for  the  f inless  version  of  the  store  shown  in  figure  6(c) 
at  incidence  in  uniform  subsonic  and  supersonic  flow.  Consider  figure  26 
showing  the  measured  normal-force  distribution  for  various  angles  of 
attack  for  Mach  nvrnber  2.  It  can  be  seen  that  considerable  loading 
exists  on  the  cylindrical  portion  of  the  store.  Also  shown  is  a solid 
line  which  is  the  force  distribution  prediction  based  on  the  slender-body 
theory  for  4°  angle  of  attack.  This  result  can  be  obtained  from  equa- 
tion (46)  of  reference  2.  The  dark  symbols  indicate  subsonic 
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experimental  data  taken  from  figure  30  in  reference  2 for  Mach  number  0.4. 
The  loading  carryover  onto  the  cylindrical  portion  appears  to  be  more 
pronounced  for  supersonic  flow  than  for  subsonic  flow.  Thus,  for  the 
store  in  subsonic  uniform  flow,  the  calculation  method  based  on  slender- 
body  theory  is  adequate.  This  is  not  the  case  for  a store  in  supersonic 
uniform  flow. 

The  pressure  distributions  along  the  windward  and  leeward  meridians 
of  the  store  are  shown  in  figure  27  as  a function  of  orifice  number  for 
5°  angle  of  attack  and  Mach  number  2.  The  following  sketch  serves  to 
relate  orifice  number  to  distance  along  the  pressure  distribution  model 
which  has  the  same  shape  and  dimensions  as  the  f inless  store  under 
consideration.  Orifice  number  6 lies  very  close  to  the  shoulder.  The 


All  dimensions  in  inches 


measurements  Hhow  loading  carryover  onto  the  cyl indrical  part  of  the 
store,  especially  on  the  leeward  side.  Theoretical  results  were  obtained 
with  a special  version  of  the  axisymmatric,  supersonic  body  modeling 
program  and  are  seen  to  overshoot  the  measured  pressures  at  the  shoulder. 
Agreement  between  theory  and  experiment  at  all  other  locations  is  good. 
Thus,  at  least  for  the  case  of  uniform  flow,  the  three-dimensional  flow 
model  presently  used  for  the  fuselage  and  store  in  supersonic  flow  is 
quite  capable  of  predicting  pressures  on  the  store  surface. 
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6.  CONCLUDING  REMARKS 


This  report  presents  the  results  of  an  effort  to  extend  the  sub- 
sonic store  separation  prediction  method  of  reference  1 to  supersonic 
speeds.  The  results  reported  herein  were  obtained  during  a combined 
theoretical/experimental  research  program,  in  order  to  isolate  and 
identify  problem  areas  in  the  flow  modeling  and  store  load  calculation 
methods,  the  configurations  of  interest  are  simple  ones  which  utilize 
wings  without  dihedral  or  incidence  attached  to  axisymmetric  fuselages. 

A pylon  can  be  attached  to  the  wing  or  fuselage  and  the  wing-fuselage 
junction  can  be  located  above  or  below  the  midplane.  Both  pylon  and 
wing  may  have  thickness.  As  one  of  the  results  of  this  effort,  a 
computer  program  for  this  class  of  configurations  was  written  with  a 
view  towards  building  it  up  to  the  general  capabilities  of  the  subsonic 
store  separation  program  of  reference  1. 

The  primary  purpose  of  this  effort  was  to  apply  linear,  potential 
flow  modeling  methods  to  the  aircraft  components  and  to  develop  methods 
for  predicting  the  flow  fields  in  which  the  store  would  be  immersed. 

The  existing  store  loading  calculation  methods  developed  in  reference  2 
would  then  be  employed.  Shortcomings  of  either  the  flow  field  or 
loading  calculation  methods  were  to  be  determined  by  comparing  with 
experimental  data.  Necessary  corrections  to  the  theoretical  methods 
were  to  be  made  to  improve  agreement,  especially  to  the  flow  field 
calculation  methods. 

The  flow  models  for  the  aircraft  components  and  the  methods  used  to 
account  for  mutual  interference  between  the  components  are  based  on 
linear,  potential  theory  for  supersonic  flow.  Line  sources,  sinks  and 
line  doublets  are  distributed  along  the  fuselage  centerline  to  represent 
volume  and  angle  of  attack  effects.  A layout  of  constant  u-veloc.ity 
panels  on  the  wing  and  pylon  (if  present)  model  lift  including  upwash 
induced  by  the  fuselage,  and  an  additional  layout  of  source  panels  model 
wing  and  pylon  thickness  effects.  Wing-pylon  mutual  interference  is 
fully  accounted  for.  Wing-to-fuaelage  interference  is  handled  by  constant 
u~velocity  panels  in  an  interference  shell  placed  over  that  part  of  the 
fuselage  influenced  by  the  wing. 

In  conjunction  with  the  analytical  work,  an  extensive  and  systematic 
wind-tunnel  test  program  was  carried  out  to  provide  flow  field  data  for 
checking  and  improving  the  accuracy  of  the  flow  models  and  the  flow  field 
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calculation  Methods.  In  addition,  force  distributions  and  overall  force 
and  moment  data  were  obtained  for  the  store  in  order  to  assess  the  adequacy 
at  supersonic  speeds  of  the  present  store  load  calculation  methods. 

After  the  computer  program  based  on  linear,  potential  flow  theory 
was  written,  it  was  necessary  to  make  corrections  to  the  method.  The 
corrections  account  for  the  fact  that  the  disturbance  from  the  wing 
leading  or  wing  root  leading  edge  (the  shock)  does  not  coincide  with  the 
Mach  lines  associated  with  linear  theory.  Also,  simple  shock  theory  was 
used  to  develop  a limit  on  the  streamwise  thickness  slope  permissible 
with  blunt  leading  edges.  Fair  to  very  good  agreement  was  obtained 
between  experiment  and  prediction  for  flow  fields.  It  was  found  necessary 
to  smooth  the  predicted  flow  field  because  the  calculated  flow  field  is 
basically  not  smooth.  The  salient  features  of  supersonic  flow  under  the 
wing  and  fuselage  were  reproduced.  With  a pylon  attached,  the  flow  is 
influenced  heavily,  and  the  present  flow  calculation  method  predicts  the 
pylon  effects  well.  The  only  systematic  deviation  occurred  in  the  region 
directly  behind  the  shock  and  it  became  more  pronounced  as  the  free-stream 
Mach  number  and  angle  of  attack  increased. 


When  the  store  is  introduced  into  the  flow  field,  additional  pertur- 
bations are  added  to  the  flow  caused  by  multiple  shock  reflections  taking 
place  between  the  store  and  the  wing.  A better  method  is  needed  to 
determine  the  effects  of  this  type  of  wing-store  interference.  The 
present  loading  calculation  method  for  the  store  body  needs  to  be  made 
three  dimensional  since  the  axial  velocity  gradients  occurring  in  the 
flow  are  basically  too  large  to  obtain  accurate  results  with  slander- 
body  theory.  This  is  in  contrast  with  subsonic  flow  where  flow  features 
are  much  more  softened  and  the  velocity  gradients  are  much  smaller  in 
magnitude.  Despite  it3  limitations  regarding  store  loading  distribu- 
tions, the  method  predicted  generally  good  results  for  total  forces  and 
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used  for  trajectory  predictions. 


In  the  following  section,  recommendations  are  described  to  improve 
the  store  load  calculation  method  and  to  further  improve  the  flow  field 
predictions. 
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7 . RECOMMENDATIONS 


The  following  recommendations  are  made  to  increase  the  accuracy  of 
the  computer  program  written  as  part  of  the  work  performed  during  the 
present  program. 

7.1  Store  Load  Calculation  Method 

The  present  method  is  capable  of  predicting  supersonic  flow  fields 
under  the  wing  or  fuselage  of  a basic  configuration.  However,  the  store 
load  distribution  calculation  method  should  be  improved.  This  can  be 
accomplished  while  at  the  same  time  improving  the  calculation  of  the 
flow  field  in  which  the  store  is  actually  immersed  by  accounting  better 
for  store  nose  shock  reflection  from  the  wing. 

A three-dimensional  approach  should  be  developed  using  the  super- 
sonic line  sources,  sinks  and  line  doublets  to  compute  the  load  distri- 
bution along  the  store  body.  Such  a method  was  developed  for  the 
subsonic  case  in  reference  2 and  the  upgraded  method  was  successful  in 
appreciably  improving  loading  results  predicted  by  the  simpler,  two- 
dimensional  loading  method  for  high  angles  of  attack.  In  contrast  with 
supersonic  flow,  the  simpler  method  was  sufficiently  accurate  in  most 
cases  encountered  for  subsonic  flow  and  was  therefore  not  implemented 
in  the  subsonic  computer  program.  The  smoothing  procedure  applied  to 
the  flow  field  in  which  the  store  is  immersed  can  probably  be  eliminated 
if  the  three-dimensional  approach  based  on  imaging  of  the  store  relative 
to  the  wing  is  used  as  described  below. 

in  addition,  the  three-dimensional  approach  should  be  accompanied 
by  an  upgraded  buoyancy  force  calculation  method.  Buoyancy  forces  can 
be  calculated  on  the  basis  of  obtaining  pressures  (using  the  Bernoulli 
pressure/velocity  relationship)  at  a number  of  points  equally  spaced  on 
the  body  circumference  in  the  crossflow  plane  without  the  store  effects 
included  in  the  flow  field.  The  buoyancy  forces  can  then  be  determined 
through  an  integration  scheme. 

7 . 2 Store-to-Aircraft  Interference 

The  next  recommendation  is  concerned  with  an  improvement  in  the 
store-to-aircraft  interference  account,  presently,  this  interference 
is  accounted  for  by  including  store  effects  in  the  flow  tangency  boundary 
condition  applied  at  the  control  points  of  the  constant  u-velocity  panels 
on  the  wing,  pylon,  and  body  interference  shell.  In  order  to  accurately 


model  the  multiple  reflections  of  the  store  nose  shock,  a very  large 
number  of  such  panels  must  be  distributed  over  the  affected  areas  of 
the  aircraft.  Furthermore,  it  would  be  necessary  to  recalculate  the 
singularity  strengths  representing  the  wing-pylon-fuselage  as  the  store 
moves  away.  In  the  improved  method  the  wing  or  the  fuselage  under  which 
the  store  is  located  is  treated  as  a reflection  plane;  that  is,  an  image 
store  is  placed  on  the  other  side  of  the  wing  or  inside  the  fuselage. 

As  the  store  drops  away,  the  singularity  strengths  representing  the 
parent  aircraft  remain  the  same  and  the  imaging  system  changes  so  that 
changes  in  the  wing-store  interference  can  be  accounted  for  without  the 
extensive  computer  time  which  would  be  required  to  recalculate  the  entire 
flow  field  at  each  step  in  the  trajectory. 

If  the  store  moves  longitudinally  from  the  wing,  then  the  imaging 
scheme  must  be  modified  because  the  wing  is  not  infinite  in  extent  and 
will  act  as  a reflection  plane  only  for  part  of  the  store. 

7.3  Further  Improvements  in  Flow  Field  Prediction  Methods 

As  an  additional  improvement  to  the  present  method,  it  is  possibxe 
to  calculate  the  flow  field  with  more  accuracy  by  replacing  the  straight 
line  characteristics  of  the  linear  theory  with  a good  approximation  of 
the  exact  (curved)  ones.  A sufficiently  good  approximation  to  t;he 
strengths  of  the  disturbances  propagated  along  the  exact  characteristics 
can  still  be  determined  from  linear  theory.  It  is  possible  to  employ 
available  theories  for  the  purpose  of  relating  flow  properties  along 
exact  characteristics  to  those  along  approximate  characteristics.  Such 
theories  are  described  in  references  18  and  19. 


NOTE:  The  transformation  in  y-axes  requires  change 

in  sign  of  sidewash  after  superposition 
is  performed. 

Trapezoidal  Panel  « (Triangle  2 - Triangle  1)  - (Triangle  4 - Triangle  3) 
(a)  Panel  with  swept  back  leading  and  trailing  edges  on  left  wing. 


neg. 


Strip  with  Swept  Back  Leading  Edge  = (Triangle  1 - Triangle  2), 

(b)  Panel  edge  with  forward  sweep  for  panel  on  left  wing. 
Figure  4.~  Superposition  scheme  for  panels  on  the  left  wing. 
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(b)  Normal- force  distribution. 
Figure  5.-  Concluded. 
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All  dimensions  in  inches 


(b)  Details  and  dimensions  of  the  pylons. 


Figure  6.-  Continued. 
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Figure  7.-  Concluded 
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Figure  8 Concluded 
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(c)  Upwash. 


Figure  10.-  Concluded. 
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Figure  11.-  Flow  field  under  the  wing  at  angle  of  attack  at 
the  one-third  semiapan  location,  a:„  - 1.37  inches. 
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Figure  12 Flow  field  under  the  wing  at  angle  of  attack  at 
the  one-third  semispan  location,  z^  - 2.12  inches. 
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(a)  Sid awash . 

Figure  13.-  Flow  field  under  the  wing  at  angle  of  attack  at 
the  one-third  samiapan  location,  zQ  - 4.37  inches. 
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Figure  13.-  Concluded 


(b)  Upwash. 
Figure  14.-  Concluded. 
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Figure  15.-  Effect  of  the  pylon  on  the  flow  field  under  the  wing 
at  angle  of  attach  at  the  one-third  semispan  location. 
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Figure  15.-  Concluded. 
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Figure  19.-  Effect  of  the  pylon  on  the  flow  field  under  the  wing 
at  the  one-third  semispan  location,  a-  « 0°. 
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Figure  22.-  Flow  field  under  the  wing  at  the 
one-third  semispan  location,  af  “ 5°. 


Figure  2 7.-  pressure  distributions  along  upper  and  lower  meridians  of  the 
pressure  distribution  store  model  in  uniform  flow. 


APPENDIX  I 

BODY  OF  REVOLUTION  IN  SUPERSONIC  FLOW 

1-1.  INTRODUCTION 

Although  the  basic  underlying  theory  is  well  known,  this  appendix 
summarizes  the  theoretical  analysis  used  to  model  the  flow  around  an  axi- 
symmetric  body  in  supersonic  flow  for  the  sake  of  completeness.  The 
method  makes  use  of  linear,  potential  flow  theory.  It  is  shown  in  the 
main  text  of  this  report  that  the  problem  can  be  split  in  two  separate 
problems:  a body  of  revolution  in  axial  flow  and  a body  of  revolution 

in  crossflow.  The  following  discussion  is  concerned  with  the  type  of 
singularities  chosen  for  each  problem  and  the  determination  of  the 
constants  in  the  distribution  of  the  strengths  of  the  singularities  from 
the  flow  tangency  boundary  conditions.  References  are  made  to  equation 
numbers  in  the  main  text  of  this  report. 

The  theory  described  here  has  been  collected  from  references  6 
and  7.  Reference  8 describes  an  early  application  of  the  axial  flow 
part  of  the  solution  to  the  calculation  of  pressures  on  the  nose  of  a 
projectile. 

1-2.  CHOICE  OF  SINGULARITIES 
1-2.1  The  Axial  Flow  Problem 

The  axial  flow  potential  is  governed  by  equation  (5)  after 

omitting  the  angular  dependence  term. 
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The  solution  for  this  potential  may  be  represented  as  an  integral  ov^r 
an  unknown  distribution  of  sources  along  the  body  axis. 
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The  unknown  function  f(£)  will  be  described  later  and  £ is  a dummy 
variable.  The  meaning  of  the  upper  limit  of  integration  is  illustrated 
in  the  sketch  below.  The  sources  are  distributed  along  the  xfi  axis 


from  0 to  L.  To  obtain  the  value  of 


at  P(V  rR) , it  is  only 


Mach  cone 
from  origin 


,(VrB> 


,Mach  cone 
from  £ 


£“xB  '3rB  -V 

necessary  to  include  sources  only  up  to  - prR  since  the  sources 

downstream  of  this  point  have  no  effect  on  the  conditions  at  P(xB,rB) . 
The  perturbation  velocities  associated  with  the  axial  potential  are 
obtained  by  taking  derivatives  of  equation  (1-2)  in  accordance  with 
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(1-4) 


Component  Ug  a i-a  the  perturbation  velocity  along  the  xfl  axis  and 

v„  is  positive  outwards  along  a radial  direction, 
u , a 

Since  the  variables  xB  and  rR  appear  in  the  upper  limit  and  the 
integrand  in  equation  (1-2)  becomes  singular  at  the  upper  limit,  it  is 
convenient  to  change  the  variable  of  integration  and  apply  Leibnitz's 
rule  for  differentiation  under  the  integral  sign.  Following  the  procedure 
described  in  reference  6,  put 
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Then  the  axial  potential  can  be  rewritten  as 
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Consequently,  from  equation  (1-4)  the  perturbation  velocities  induced 
by  the  axial  potential  are  given  by 
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where  f'  denotes  differentiation  of  function  f with  respect  to  its 
argument.  It  has  been  assumed  that  f (0)  = 0.  Later  it  will  be  shown 
that  this  condition  corresponds  to  bodies  with  pointed  noses. 

A particular  simple  and  useful  function  is  the  one  associated  with 
a line  source. 


f(£)  **  K£ 
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Substituting  this  function  into  equation  (1-6)  and  using  equation  (1-5) 
for  £ gives 
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The  associated  perturbation  velocities  can  be  found  by  direct  differen- 
tiation in  accordance  with  equation  (1-4) . 
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Both  nonzero  velocities  are  functions  of  parameter  x_/6r-,.  Thus, 

they  are  invariant  along  lines  for  which  xg/rB  is  constant;  that  is, 

along  rays  from  the  origin.  The  flow  field  described  by  equation  (1-10) 

is  called  conical.  For  a cone  with  vertex  angle  6 . reference  6 

cone' 

relates  the  constant  K to  the  angle  through 


reference  6 
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The  conical  flow  field  is  shown  in  the  following  sketch.  Velocities 
Ug  a and  vB  a are  constant  along  the  rays  from  the  origin.  On  the  Mach 
cone  the  perturbation  velocities  are  zero  and  the  flow  is  still  parallel 
to  the  free  stream.  Two  streamlines  arc  shown,  one  originating  in  the 
free  stream  and  one  emanating  from  a point  on  the  line  source  along  the 
axis.  Any  streamline  can  represent  a solid  surface,  but  there  is  only 


Mach  cone 


Streamlines 


Origin 

one  streamline  coincident  with  a ray.  A cone  containing  that  ray  repre- 
sents a solid  cone  in  supersonic  flow  at  least  within  the  constraints  of 
linear  theory.  The  flow  field  associated  witn  this  cone  is  specified 
by  equation  (1-10)  with  equation  (1-11) . 

The  axial  flow  around  an  arbitrary  but  pointed  body  of  revolution 
aligned  with  the  flow  can  be  obtained  by  a superposition  of  conical 
solutions  of  the  type  discussed  above  since  the  starting  procedure  is 
clear.  The  procedure  is  indicated  in  the  following  sketch.  For  the 
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ogive  shape  shown  below,  the  composite  source  distribution  for  the 
axial  flow  solution  is  obtained  by  generating  a superposition  of  linearly 
varying  sources;  that  is,  a superposition  of  cone  solutions. 

f(0  - k0$  - ex)  +k2(^  - is)  + ...  d-12) 


Constants  K^Kg,  and  K are  negative  for  the  example  in  the  sketch  above, 

and  £Q  “ “ XS1  ~ PrB1  > ia  “ xBa  “ PrB£>  etc,,  are  the  origins  of 

the  line  singularities,  The  constants  can  be  evaluated  in  a step-by-step 

manner  as  follows.  The  constant  Kq  corresponds  to  the  solution  which 

satisfies  the  flow  tangency  condition  at  (xB  ,rB  ) and  gives  the  flow  near 

the  nose  tip.  At  the  next  control  point  on  the  body,  it  is  necessary  to 

superimpose  another  conical  solution  with  its  origin  located  at 

£ ■ xn  - pr„  and  its  constant  K is  chosen  so  that  the  flow  tangency 

i “1  Bi  1 

condition  on  the  body  at  point  (xg^rg^)  is  satisfied.  The  second  solution 

does  not  affect  the  nose  flow  ahead  of  the  Mach  cone  with  its  origin  at 

£ . It  should  be  noted  here  that  this  particular  property  is  charaeter- 
i 

istic  of  solutions  to  the  wave  equation  given  by  (i-l) . Let  NXBOPY  be  the 
number  of  line  sources  with  origins  equally  spaced  on  the  body  centerline. 
As  a result,  the  number  of  control  points  equals  NXBODY. 
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The  superposition  scheme  for  representing  a body  of  revolution  in 
supersonic  flow  is  applicable  only  if  the  body  meridianal  section  is 
smooth.  For  bodies  with  corners,  additional  solutions,  such  as  Rrandtl- 
Meyer  expansion  theory,  are  required  to  properly  account  for  the  flow 
near  the  shoulders.  The  layout  of  the  body  definition  points  and  control 
points  as  shown  in  the  proceeding  sketch  is  used  in  the  body  flow  modeling 
program  of  references  9 and  20.  Following  a suggestion  by  R.  C. 

Carmichael  of  NASA/Ames  Research  Center  associated  with  reference  21, 
the  control  points  can  be  shifted  forward  a half-space,  AxQ/2 . This  modi- 
fication tends  to  smooth  the  solution.  The  modified  layout  of  the  body 
definition  and  control  points  is  indicated  in  the  following  sketch.  In 
this  case,  the  number  of  line  sources  NXBODY  equals  the  number  of  control 
points . 


The  flow  tangency  boundary  condition  for  the  axial  flow  problem  given 
by  equation  (8a)  can  be  rewritten  as  follows. 


dK 

dxB 


- x 


first  control 
point 


U-13) 


Substituting  from  equation  (1-10) 


then  gives 
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XB  " XB1  “ (Axb/2) 
rB  - rBi  - R(xBi) 

where  x " and  rn  are  the  coordinates  of  the  first  control  point. 

B1 

With  Kq  known,  the  second  constant  can  now  he  explicitly  determined 

from  boundary  condition  (8a)  applied  to  the  second  control  point  (xB  ,rR  ) 

2 £ 

on  the  body  surface. 


(1-15) 


The  terms  in  the  above  expression  are  given  by 


£ - 0 
so 


^ " xBd1  “ ^Bd, 


xB  - xB  + 


Ax  „ body  length 
B : XBODY 


(1-16) 


(Continued  on  next  page) 


• « 


specified 


The  function  f'  ({•)  is,  go  far,  unknown  and  is  to  be  determined  from  the 
flow  tangency  condition  in  the  crossflow  plane  given  by  equation  (8b) . 
To  avoid  confusion  with  the  function  f(£)  used  in  the  social  solution, 


d(£)  - f ($) 


(1-21) 


since  we  are  dealing  now  with  a doublet  distribution. 

In  accordance  with  the  procedure  shown  in  reference  7 , the 
following  linearly  varying  distribution  is  used  for  the  doublet  strength 
distribution  as  given  by 


d(£)  - Kd£ 


(1-22) 


Again,  letting  £ " ^ " P^coah  a as  in  equation  (1-5),  and  carrying 
out  the  integration  indicated  by  (1-20)  results  in  a three-dimensional 
doublet  solution 


VW3*  “ Kdcos  6 


2rB  V 


*B  ~ 3 


zra 

T. 


(1-23) 


Direct  differentiation  of  the  crossflow  potential  gives  the  pertur- 
bation velocities  associated  with  the  line  doublet. 


Vd  V 

V “ VW—  “ cos  0 

v»  oXjj  d v p2r2 


- 1 
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B^d  V 
SrT 

B 


V, 


P2cos  0 


oo»h-1  ( ds-)  + i-  - 1 

0rB  V02^ 


(1-24) 


Vd  . 1 


K.sin  0 
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The  crossflow  around  an  arbitrary  body  of  revolution  can  be  obtained 
by  a superposition  of  doublet  solutions  in  the  same  manner  used  for  the 
axial  flow  problem.  Thus,  the  composite  doublet  distribution  for  the 
crossflow  solution  is  obtained  by  generating  a superposition  of  linearly 
varying  doublet  distributions 


d<4) 


<d,o^  +Vr^  “ V + 


Kd,e(?  - 


,)  + 


(1-25) 
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in  terms  of  constants  Kd , q,k6> 1 ,Kd> a and  £0  - 0,  ^ - xBdi  - arBdi , 

■ xBd2  ~ ^rBd2i  etc.,  are  the  origins  of  the  line  singularities.  These 
constants  are  also  evaluated  in  the  same  step-by-step  manner  employed  in 
the  axial  flow  problem  but  with  the  boundary  condition  specified  by 
equation  (8b) . It  is  noted  that  the  perturbation  velocity  terms  required 
by  that  equation  are  given  by  the  first  two  of  equation  (1-24)  and  contain 
the  term  cos  0.  Consequently,  the  cos  8 term  is  calcelled  in  equa- 
tion (8b) . The  first  constant  is  obtained  from  the  boundary  condition 
applied  at  the  first  control  point  Rewriting  equation  (8b) 

gives 


a^cos  9 


- **  1 


,d  dR  qj  ,d 

Voo  d*B  " V„ 


Substitution  from  equation  (1-24)  leads  to 


(1-26) 


Ho' 


^2  cosh 


1 6‘ 
J i 


(1-27) 


Let  NXBODY  be  the  number  of  line  doublets  with  origins  equally  spaced  on 
the  body  centerline.  In  general,  the  ntk  constant  is  explicitly  d ;ter- 
mined  from  the  boundary  condition  applied  at  the  (n+l)th  control  point 
as  follows. 


af  " Kd, k-i 


f>  rc 03h-x  i ~ M + (\^  - ^ 

) 2 \ PrB  . V PrB+  J 

\ l.  v n+i  ' \ n+i  / 


"Bn+1  " 5k-l} 


■ 

dxB 


n+i  


(1-28) 


(Continued  on  next  page) 
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d,n  b 


where 


n - 1,2, ... (NXBODY-1) 


^o  * 0 


$k-i  ’ XBd 


k-x  * " " 2*3> 


(Concluded) 


xBd  ” ^rBd 

n n 

xR  - x + ax 

“n+i  Bn  S (1-29) 

. „ Body  length 

AH3  “ NXBODY 

rR  “ R(xh  ) 

Bn+i  Sn+i 

dir  - r,<xb  5 

B n+x  n+1  J 

keferences  9 and  20  give  an  equivalent  form  of  equation  (1-28) . 

T >3.  BODY  INDUCED  PERTURBATION  VELOCITIES 

The  induced  perturbation  velocities  due  to  one  line  source  are  given 
by  equation  (1-10) . The  contributions  from  one  line  doublet  are  specified 
by  equation  (1-24) . The  perturbation  velocities  induced  by  a body  with 
a number,  NXBODY,  of  line  sources  and  doublets  distributed  along  the 
centerline  are  obtained  by  adding  the  effects  of  all  the  singularities. 
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The  constants  K and  K,  have  been  determined  above.  At  a field  point 
n a,n 

specified  in  body  coordinates  xB>rB>®>  the  body  induced  perturbation 
velocities  are 


xBd  *rBd  " coordinates  of  the  n body  definition  point 
n n 

Ug  * axial  perturbation  velocity  component 

v_  - radial  perturbation  velocity  component 
& 

wQ  « tangential  perturbation  velocity  component 

The  positive  directions  are  shown  again  for  convenience  in  the  following 
sketch  - 


APPENDIX  II 


I 

| PERTURBATION  VELOCITIES  INDUCED  BY  A 

[ SEMI-INFINITE  TRIANGLE 


I I - 1 . INTRODUCTION 

This  appendix  describes  the  perturbation  velocities  induced  by  a 
j semi-infinite  triangle  with  a sweptback  leading  edge,  as  shown  in  figure  3 

of  the  main  text  in  supersonic  flow  for  two  cases.  First,  the  triangle 
is  subjected  to  a constant  jump  in  axial  velocity  across  its  plane. 

Second,  a source  distribution  of  constant  strength  is  placed  in  its  plane. 
Even  though  certain  elements  and  difficulties  of  the  underlying  theory  are 
discussed,  the  primary  purpose  of  this  appendix  is  to  specify  the  formula- 
tions for  perturbation  velocities  as  used  by  the  computer  program. 

II-2 . ANALYTICAL  APPROACH 

The  general  method  of  finding  the  potential  flow  solution  for  a 
semi-infinite  triangle  with  a swept  back  leading  edge  in  supersonic  flow 
is  given  in  some  detail  in  references  9 and  20.  A summary  of  the  method 
is  contained  in  reference  10.  In  particular,  the  last  reference  discusses 
the  difficulties  associated  with  the  nonplanar  case.  For  that  case,  ■'he 
plane  of  the  semi-infinite  triangle  makes  an  angle  with  the  x-axis  in 
contradistinction  with  the  planar  case  shown  in  figure  3.  It  is  shown 
in  reference  10  that  the  solution  obtained  in  references  9 or  20  for 
the  nonplanar  case  can  introduce  additional  streamwise  vortex  sheets  in  the 
flow  field.  In  addition,  the  nonplanar  solution  features  a built-in  thick- 
ness distribution  for  the  case  of  the  triangle  being  subjected  to  a jump 
condition  in  u-velocity.  In  principle,  this  undesirable  thickness  distri- 
bution is  removable  by  subtracting  another  thickness  distribution  to  cancel 
the  effects  of  the  undesired  one.  However,  a simple  means  for  eliminating 
the  unwanted  vortex  sheet  of  vorticitv  has  not  been  found.  Accordingly, 
only  planar  solutions  are  used  for  the  case  of  triangles  with  a jump  in 
axial  velocity.  For  simplicity,  only  planar  solutions  for  triangles  with 
constant  source  strength  are  employed  herein. 

For  the  sake  of  completeness,  the  method  of  solution  will  now  be 
outlined.  The  linearized  differential  equation  for  the  velocity  potential 
generated  by  a small  perturbation  of  a steady  supersonic  flow  characterized 
by  Mach  number,  ^ , is  given  by  the  three-dimensional  wave  equation 
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Perturbation  velocities  u,v,w  also  obey  this  wave  equation.  The  coordi- 
nate system  is  shown  in  figure  3.  Let  12  be  an  arbitrary  variable  repre- 
senting either  the  potential  0 or  any  one  of  the  perturbation  velocities. 
Reference  22  describes  how  a general  solution  to  equation  (II-l)  can  be 
obtained  using  Volterra's  solution  to  the  two-dimensional  version  of  the 
wave  equation.  Considering  a small  flow  perturbation  originating  on  a 
surface,  S,  the  value  of  either  the  potential  or  a perturbation  velocity 
at  any  point,  P(x,y,z),  is  given  by  the  following  integral  equation. 


S2(x,y  ,z) 


1 S 


2ir  "5x 


<y  dS 


+ 3¥S  //  ® - «')  tds 

T 


(11-2) 


The  domain  of  dependence  x is  the  Mach  forecone  from  point  P.  Both 
integrals  are  to  be  evaluated  over  the  part  of  surface  S intercepted  by 
the  Mach  forecone  from  P.  The  unprimed  variable  denotes  the  value  of 
that  variable  on  the  same  side  of  surface  S as  point  P while  the 
primed  variable  denotes  its  value  on  the  opposite  side  of  S.  The  function 
0 is  the  particular  solution  of  equation  (II-l)  specified  by  Volterra 
which  vanishes,  together  with  its  derivative  dcr/dv , everywhere  on  the  Mach 
forecone  from  P. 


a - cosh-1  — - ~ ^ — (11-3) 

pV(y  - 7)) 2 + (z  - 02 

In  fact,  this  function  is  the  indefinite  integral  of  the  fundamental  solu- 
tion to  the  wave  equation  (II-l)  associated  with  a three-dimensional  source 
at  (£>n>C)  i-n  supersonic  flow.  The  variable  v is  associated  with  a 
vector  called  the  conormal  to  the  surface  S.  This  vector  is  related  to 
I the  normal  vector  to  that  surface.  It  has  a special  property  in  that  the 

conormal  to  the  Mach  forecone  lies  in  the  surface  of  the  Mach  cone. 
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After  performing  substitutions  for  a.da/dv,  and  d£  in  equation  (II-l), 
it  is  possible  to  rewrite  the  integrals  in  equation  (II-2)  for  either  an 
inclined  triangle  or  a triangle  in  the  x-y  plane  such  as  the  one  shown 
in  figure  3.  The  integrals  can  then  be  evaluated  if  the  terms  1,0  - O') 
and  (dft/dv  + dft'/dv')  are  prescribed  on  surface  s which  in  this  case  is 
the  semi-infinite  triangl6.  It  is  most  convenient  (but  possibly  somewhat 
restrictive)  to  set  the  terms  equal  to  a constant  or  aero.  One  choice  is 
to  specify  a constant  jump  in  axial  velocity  u across  the  triangular 
surface  corresponding  to  a lifting  surface.  Another  choice  is  to  set  12 
equal  to  the  potential  4>  and  the  result  is  the  integral  form  of  the  poten- 
tial due  to  a surface  distribution  of  sources. 

In  the  following  discussion,  the  formulations  have  been  especially 
tailored  to  suit  the  computer  program.  They  have  been  deduced  from  the 
results  obtained  by  Woodward  as  published  in  reference  9. 

II-2.1  Velocities  Induced  by  a Semi-Infinite  Triangle 
with  a Jump  in  u-Velocity  Across  Its  Plane 

For  this  case,  variable  12  in  equation  (II- 2)  is  set  equal  to  the 
perturbation  velocity  u+  on  the  upper  surface  of  the  semi-infinite 
triangle  shown  in  figure  3.  The  desired  solution  must  have  a constant 
discontinuity  in  u everywhere  on  S.  Reference  9 shows  the  procedure 
required  to  obtain  the  solution  expressed  in  terms  of  a pressure  jump.  In 
that  reference,  the  pressure  jump  is  actually  directly  related  to  the  axial 
velocity  jump  through  the  linear  velocity-pressure  relationship.  In  this 
work,  the  more  fundamental  constant  u-vclacity  jump  concept  is  retained. 

For  tile  reasons  mentioned  earlier  in  this  appendix,  only  the  solution  for 
the  triangle  in  the  x-y  plane  is  considered. 

The  coordinate  system  associated  with  the  semi-infinite,  sweptback 
triangle  located  in  the  x-y  plane  as  shown  in  figure  3 is  used  in  the 
following  expressions  for  the  perturbation  velocities.  In  terms  of  functions 
to  be  specified,  the  perturbati.on  velocities  at  a point  P(x,y,z)  are 
given  below  for  x > 0,  Unswept  triangles  are  discussed  separately  and 
expressions  for  the  perturbation  velocities  are  given  later. 
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The  quantity  m^e  is  the  tangent  of  the  sweep  angle,  For  x <£  0, 

the  velocities  are  zero.  The  strength  of  the  singularity  in  the  semi- 
infinite triangle  is  given  by  u^AT^.  Functions  • F4  > Fg  , and  F? 

depend  on  the  geometry  and  are  different  for  subsonie,  sonic,  and  supersoni 
leading  edges  associated  with  the  semi-infinite  triangle.  The  unswept 
leading  edge  is  dealt  with  separately.  All  four  cases  are  discussed  below 

Subsonic  Leading  Edge  m|g  > p2  , m^e  - tan  ip£e 


tan (90°-  u)  = cot  p - - 1 = p2 


z ,w 
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The  above  functions  arc  to  be  substituted  into  equation  (II-4) . On  the 
x-axis  or  along  the  leading  edge,  the  resulting  velocities  are  singular  or 
indeterminate.  The  computer  program  has  been  arranged  to  set  the  veloci- 
ties equal  to  zero  when  the  field  point  lies  either  on  the  x-axis  or  the 
leading  edge.  If  the  coordinates  of  point  P satisfy  the  relation 

x2  - (32(y2  + zE)  £ 0 (II-6) 

then  the  point  P lies  on  or  outside  the  Mach  cone  from  the  origin.  For 
this  case,  the  velocities  are  also  zero. 
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All  functions  are  the  3aine  as  for  the  subsonic  leading-edge  case  given  by 
equation  (H-5)  except  function  F^.  It  ia  now  specified  as  follows. 


- Q2  (vz  + z‘ 
x - 3y 


(I  I- 7) 


The  velocity  components  are  singular  or  indeterminate  at  points  along  the 
x-axis  or  the  leading  edge  and  are  set  equal  to  zero  by  the  program.  If 
the  coordinates  of  point  P satisfy  equation  ( XI— 6 ) , all  velocity  compo- 
nents are  zero. 


Supersonic  Leading  Edge 


mJe  < 


» mie 


- tan  f 


£e 


tan  (90  - ix)  » P* 


The  field  point  P(x,y,z)  lies  inside  the  Mach  cone  from  the  origin  when 
the  following  condition  is  satisfied  by  its  coordinates. 

x2  - P2(y2  + z2)  >0  (II-8) 

For  this  case,  all  functions  are  given  by  equation  (11-5}  except  F£ 
which  is  changed  to 


m "V  (P2  - ni2  )[x2  - p2(v2  + 

F,  *2 tan  1 ^ (ir-9) 

VP“  - nje  ^ie  - 0 y 

The  velocity  components  are  singular  or  indeterminate  at  points  along  the 
x-axis.  The  computer  program  is  arranged  to  set  u,v,w  equal  to  zero 
when  both  y and  z are  zero.  If  the  point  P(x,y,z)  lies  outside  the 
Mach  cone  from  the  origin  (equation  (II-8)  not  satisfied)  but  is  inside 
the  Mach  cone  with  its  origin  at  (x^e,v,0)  on  the  leading  edge,  then  the 
following  conditions  are  met. 
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For  this  situation,  the  functions  assume  the  following  constant  values 


If  the  point  P is  outside  the  Mach  cone  from  the  origin  and  its  y-coordi- 
nate  is  negative,  then  all  velocities  are  set  equal  to  zero.  Zero  veloc- 
ities also  result  when  the  point  is  outside  the  Mach  cone  from  the  origin 
and  the  Mach  cone  from  the  leading  edge. 

Unawept  Leading  Edge  m^e  " 0,  (32  - - 1 
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This  case  is  really  a special  version  of  the  seini-infinite  triangle  with  a 
supersonic  leading  edge.  Equation  (II-4)  does  not  apply  to  this  case. 
Instead,  the  following  expressions  specify  the  perturbation  velocities 
induced  at  a point  J?(x,y,z)  by  the  unswept  semi-infinite  triangle  with 
a constant  jump  in  u-velocity. 


J 


Functions  , F.,,  F4  and  Fy  depend  on  the  location  of  the  field  point. 

If  the  point  lies  inside  the  Mach  cone  from  the  origin,  its  coordinates 
satisfy  equation  (II-8) . Functions  F4  and  F7  are  given  by  equation  (11-5) 
and 


F - tan 

i 


-i  zVxa  - fi~(ya  + zg) 
-xy 


"s 


F0  - -[5  tan 


-l  Vx*  - fl2(y2  ± za) 


-(iy 


(11-13) 


The  velocity  components  are  singular  or  indeterminate  at  points  along  the 
x-axis.  The  computer  program  is  arranged  to  set  u,v,w  equal  to  zero 
when  both  y and  z are  zero.  If  the  point  P(x,y,z)  lies  outside  the 
Mach  cone  from  the  origin  but  is  inside  the  Mach  cone  with  its  origin  at 
(0  ,y  ,0)  on  the  leading  edge,  then  conditions  (II-6)  and  (11-10)  are 
satisfied.  The  functions  then  assume  the  following  constant  values. 
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Figure  II-l  shows  the  behavior  of  velocity  components  u/Vw,  v/Vk, 
and  w/Vm  for  unit  strength  u^/V^  along  a traverse  parallel  to  the 
y-axis  in  the  z ■ 0 plane.  Results  are  shown  for  the  subsonic,  sonic, 
and  supersonic  leading-edge  cases  and  the  Mach  number  equals  VI 

11-2.2  Velocities  Induced  by  a Semi-Infinite  Triangle  with 
a Source  Distribution  of  Constant  Strength 

For  the  purpose  of  generating  the  solution  for  this  case,  variable 
il  in  equation  (II-2)  ia  set  equal  to  the  potential  0 on  the  upper 
surface  of  the  semi-infinite  triangle  shown  in  figure  3.  The  partial 
derivative  d0/dv  in  equation  (XX-2)  represents  the  velocity  component 
along  the  conorrnal  to  the  upper  surface  of  the  triangle  and  d0‘/Dv'  is 
the  velocity  component  along  the  conormal  to  the  lower  surface,  if  0‘ 
has  the  same  sign  as  0,  a discontinuity  in  that  velocity  component  will 
appear  in  the  flow  at  the  surface  of  the  triangle.  For  the  sweptback 
triangle  located  in  the  x-y  plane  as  shown  in  figure  3 (disregarding  the 
u+  and  u_  vectors  shown) , the  velocity  component  under  consideration  is  in 
fact  the  velocity  component  w in  the  z-direction.  Therefore,  the  resulting 
solution  is  one  which  produces  an  upwash  on  the  upper  surface  of  the  triangle 
and  a downwash  just  below  the  surface.  Since  the  planar  solution  for  the 
constant  u-velocity  jump  triangle  is  used  in  the  computer  program  an 
described  in  the  preceding  section,  the  planar  solution  for  the  source 
triangle  given  in  reference  9 is  used  here  for  the  sake  of  consistency. 

The  coordinate  system  shown  in  figure  3 is  used  again  in  the  expres- 
sions for  the  perturbation  velocities  that  follow.  For  x > 0 and  with 
nonzero  leading-edge  sweep,  the  perturbation  velocities  induced  at  a given 
point  P(x,y,z)  by  a semi-infinite  triangle  with  a source  distribution  of 
constant  strength  are  given  below.  For  x <;  0 , the  velocities  are  zero. 

The  case  for  an  unswept  triangle  is  discussed  separately  and  expressions 
for  perutrbation  velocities  are  given. 
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The  source  strength,  tan  0,  of  the  semi-infinite  triangle  under  consider- 
ation is  associated  with  the  slope  of  the  thickness  envelope  of  the  surface 
under  consideration.  As  discussed  in  the  flow  field  calculation  section 
in  the  main  text,  the  thickness  slope  i3  bounded  in  magnitude.  Referring 
to  the  general  panel  layout  shown  in  figure  1,  the  required  source  strengths 
for  the  four  semi-infinite  triangles  used  to  generate  the  solution  for  one 
source  panel  shown  on  the  wing  are  given  by  the  slope  of  the  wing  thick- 
ness envelope  calculated  at  the  panel  centroid.  Functions  F1 , F£,  and  Ft 
are  determined  as  described  in  section  II-2.1  for  the  cases  of  subsonic, 
sonic,  and  supersonic  leading  edges.  In  all  three  cases,  the  v and  w 
components  are  singular  or  indeterminate  at  points  on  the  x-axis.  Further- 
more, the  u,v,and  w components  are  singular  or  indeterminate  at  the 
leading  edge.  The  computer  program  is  arranged  to  set  the  velocity  com- 
ponents equal  to  zero. 

As  already  mentioned  in  connection  with  the  constant  u-velocity  jump 
triangle,  the  unswept  leading-edge  case  is  a special  version  of  the  semi- 
infinite triangle  with  a supersonic  leading  edge.  Equation  (11-15)  is 
not  used  for  this  case.  The  expressions  below  specify  the  perturbation 
velocities  induced  at  a point  P(x,y,z)  by  the  unswept  semi-infinite 
source  triangle  for  x > 0.  If  x <[  0,  the  velocities  are  zero. 
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Functions  Fx , , and  F depend  on  the  location  of  the  point  P.  If  the 
point  lies  inside  the  Mach  cone  from  the  origin,  its  coordinates  satisfy 
equation  (II-8) . Function  F^  is  given  by  equation  (11-13)  , function  F 
by  equation  (II-5)  , but  is  changed  to 
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The  velocity  components  are  singular  or  indeterminate  at  points  along  the 
x-axis.  The  computer  program  is  arranged  to  set  the  velocities  equal  to 
zero.  It  the  point  P(x,y,z)  lies  outside  the  Mach  cone  from  the  origin 
but  is  inside  the  Mach  cone  with  its  origin  at  (0,y,0)  on  the  leading 
edge,  then  equations  (II-6)  and  (11-10)  are  satisfied.  The  functions  then 
acquire  the  following  constant  values. 


(11-18) 


It  is  seen  that  for  points  in  the  plane  of  the  semi-infinite  triangle 
between  the  Mach  line  from  the  origin  and  the  unswept  leading  edge,  side- 
wash  component  v is  zero  and  components  u and  w are  constant. 
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The  behavior  of  the  velocity  components,  u/VM,  v/V^,  and  w /v  , for 
given  source  strength  tan  0 is  shown  in  figure  II-2  for  a swept  semi- 
infinite  triangle  along  a traverse  parallel  to  the  y-axis  in  the  z «■  0 
plane.  Results  are  shown  for  the  subsonic,  sonic  and  supersonic  leading- 
edge  cases  and  the  Mach  number  is  VT 
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APPENDIX  III 
DETACHED  SHOCK  ANALYSIS 


III-l.  INTRODUCTION 

In  this  appendix,  a method  is  described  for  predicting  the  shape  of 
a detached  shock  which  stands  ahead  of  a two-dimensional  shape  with  a 
blunt  leading  edge.  For  shapes  with  cylindrical  forebodies,  Love  in 
reference  14  presents  a shock  shape  and  location  determination  method 
based  on  simple  concepts  and  augmented  with  certain  corrections  obtained 
from  experimental  data.  Recently,  Rizzi  at  the  NASA/Amea  Research  Center 
performed  calculations  to  determine  shock  standoff  distance  and  other 
properties  of  the  flow  over  a two-dimensional  flat  body  with  cylindrical 
leading  edge.  His  calculations  involve  solving  the  Euler  aquations 
(inviscid  Navier-stokes  equations)  by  means  of  a finite-differencing 
scheme  and  this  method  is  cited  in  reference  15. 

During  the  course  of  work  performed  on  the  development  of  the  flow 
field  calculation  method  described  in  this  report,  the  above-mentioned 
shock  prediction  methods  were  investigated  in  connection  with  correcting 
the  linear  flow  theory  for  nonlinear  shock  effects.  As  it  turned  out, 
the  present  flow  field  calculation  method  itself  predicts  the  shock 
location.  Nevertheless,  an  independent  shock  prediction  method  was  found 
to  be  useful  in  checking  the  shock  shape  predicted  by  the  present  method 
as  implemented  in  the  computer  program. 

In  what  follows,  the  actual  mechanics  involved  in  applying  the 
detached  shock  prediction  method  to  a cylindrical  forebody  shape  are 
described.  The  method  for  determining  the  shock  shape  is  basically 
that  of  Love  in  reference  14  but  with  standoff  distance  taken  from 
calculations  performed  by  Rizzi  of  reference  15. 

In  addition,  a criterion  for  the  maximum  value  of  the  wing  thickness 
slope  is  developed  on  the  basis  of  the  detached  shock  model.  An  upper 
bound  is  required  by  the  present  flow  field  calculation  method  when  the 
configuration  of  interest  includes  a wing  with  a blunt  leading  edge.  For 
such  a wing,  the  wing  thickness  slope  is  infinite  at  the  leading  edge. 

As  described  in  this  report,  the  strengths  of  the  source  panels  used  to 
model  thickness  are  directly  related  to  the  thickness  slopes  measured  in 
the  streamwise  direction  at  the  panel  centroids.  If  the  number  of  source 
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panels  in  a chordwise  row  is  large,  the  thickness  slopes  associated  with 
the  forward  source  panels  may  become  too  large  in  magnitude  and  are  then 
incompatible  with  the  linear,  potential  theory  on  which  the  solution  for 
the  source  panel  is  based. 

III-2.  PROCEDURE  FOR  CALCULATING  SHAPE  OF  DETACHED  SHOCK 

The  procedure  described  here  is  the  one  developed  by  Love  in  refer- 
ence 14  with  the  exception  of  the  determination  of  the  shock  detachment 
distance.  Rizzi  of  reference  15  solved  the  Euler  equations  using  a 
finite  element  technique  to  obtain  details  of  the  flow  over  a two-dimen- 
sional flat  shape  with  a cylindrical  forebody  at  transonic  or  supersonic 
speeds.  The  flow  region  was  divided  into  a mesh  of  many  quadrilateral 
panels  of  varying  size.  The  results  of  Rizzi'e  computations  in  terms  of 
shock  standoff  distance  and  sonic  line  location  are  useful  as  the  starting 
point  of  the  approximate  shock  shape  model  developed  by  Love. 

Figure  III-1  shows  the  two-dimensional  standoff  distance  of  shocks 
associated  with  flow  normal  to  a circular  cylinder  as  a function  of  Mach 
number,  M^ , normal  to  the  body.  The  theoretical  points  represent  the  con- 
verged results  from  calculations  performed  by  Rizzi  for  different  mesh 
sizes.  Also  shown  are  a few  points  from  experimental  measurements  published 
in  reference  6.  Agreement  between  theory  and  experiment  ia  very  good 
for  the  two  experimental  points  shown.  Rizzi's  calculations  also  include 
the  location  of  the  sonic  point  on  the  body.  This  location  is  schemati- 
cally shown  in  figure  III-2.  Although  the  calculated  standoff  distances 
showed  convergence  with  finer  mesh  size  used  by  Rizzi,  the  sonic  point 
locations  did  not.  Apparently,  the  angular  locations  of  the  sonic  point 
require  even  finer  or  differently  spaced  quadrilaterals  before  convergence 
is  reached.  Work  is  presently  under  way  to  try  a meBh  size  which  is 
exponentially  varied  and  is  very  fine  near  the  body  surface. 

The  equations  required  for  calculating  the  shock  shape  by  the  approx- 
imate method  of  Love  are  given  in  reference  14  and  will  now  be  expressed 
in  the  coordinates  (x,y)  and  radius  d/2  shown  in  figure  iTI-2.  Based  on 
a hyperbolic  representation  of  the  shock,  the  shape  of  the  shock  is  given 
as  follows. 
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The  ratio  xQ/d'  is  given  by  the  following  expression: 


(3atan  e - p'\/^atan£e  - 1 + tan  tj 
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Angle  is  the  shock  wave  angle  for  sonic  flow  behind  the  shock.  It. 

can  be  calculated  from  equation  (167)  in  reference  11. 


(III-3) 


Angle  ij  is  the  angle  between  the  normal  to  the  free  stream  and  the 
control  line.  The  most  convenient  and  accurate  way  to  determine  it  is  to 
read  its  value  from  figure  10  in  reference  14  as  a function  of  M • 
Finally,  quantity  x’/d'  is  related  to  the  shock  standoff  distance 
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In  this  equation,  AxR/d  is  given  in  figure  Ill-i  as  a function  of 
Angle  6det  is  the  semiapex  angle  for  a wedge  for  which  the  shock  is 
first  detached.  It  can  be  determined  from  chart  2 in  reference  11, 
or  the  following  relations  can  be  used.  The  shock  wave  angle  f°r 

maximum  stream  deflection  behind  the  shock  (which  occurs  when  it  is  first 
detached)  is  given  by  equation  (168)  in  reference  11. 
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The  corresponding  maximum  wedge  angle,  <5det,  or  the  angle  for  which  the 
shock  first  is  detached  is  given  by  equation  (138)  in  reference  11. 
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I I 1-3.  MAXIMUM  ELOPE  ANGLE 

The  actual  flow  behind  the  detached  shock  includes  the  following 
features.  A region  of  subsonic  flow  exists  between  the  body  surface,  the 
shock,  and  the  sonic  line  which  is  curved  in  reality.  At  the  sonic  line, 
the  flow  has  accelerated  from  subsonic  and  enters  the  supersonic  speed 
region.  Ho  part  downstream  of  the  sonic  point  on  the  body  affects  the 
detachment  process.  The  situation  is  depicted  schematically  in  the 
following  sketch. 


From  the  sonic  point  on,  the  flow  accelerates  to  uniform  speed  in  the 
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lines  (characteristics)  originating  at  points  on  the  body  behind  the  sonic 
point  are  propagated  out  to  the  shock  and  causa  it  to  curve. 


On  the  basis  of  the  above  observations,  the  first  body  slope  used  in 
the  supersonic  flow  field  calculation  method  is  bounded  by  the  slope  of 
the  body  at  the  sonic  point.  Until  Rizzi's  calculations  generate  converged 
results  for  the  values  of  angle  @B  associated  with  the  sonic  point,  angle 
^det  should  be  used  as  the  upper  limit  on  the  surface  slopes.  This  result 
is  based  on  the  assumption  that  the  nose  contour  in  front  of  the  sonic  line 
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will  not  change  the  detached  shock  location.  This  assumption  is  satis- 
fied if  the  nose  contour  lies  within  the  wedge  formed  by  drawing  the 
tangent  to  the  nose  at  the  sonic  point.  This  assumption  was  first 
suggested  by  Busemann  (ref.  23) . As  can  be  seen  from  figure  I1I-3, 
angle  6det  approaches  aero  as  the  Mach  number  is  reduced  to  1. 

In  the  application  of  the  above  to  blunt  wings  with  marginally 
supersonic  leading  edges,  it  is  probably  better  to  consider  the  edge 
subsonic.  The  upper  limit  on  the  thickness  slopes  is  then  determined 
on  the  basis  of  the  free-stream  Mach  number  as  described  in  this  report. 


Shock  standoff  distance  for  cylindrical  forebodies. 


Figure  III-2.-  Coordinates  associated  with  detached  shock  model 


Figure  III-3.-  Variation  of  shock  detachment  angle  with 
upstream  Mach  number  for  wedges. 
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